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specific  vibrational  levels  of  the  electronic  states  corresponding  to  the 
strongest  radiating  band  systems  belonging  to  the  uranium/oxygen  and 
titanium/oxygen  systems  and  prediction  of  TR  and  possible  optical  oscillator 
st  rengths . 


The  results  obtained  to  date  for  the  uranium/oxygen  system  indicate  that 
both  U0'K  and  UO  ,+  are  very  strong  LWIR  radiators  in  the  10-14  inn  region.  In 
addition,  U0+ ,  which  is  a  potential  late  time  radiator,  exhibits  strong 
visible  bands  and  should  be  efficiently  solar  pumped.  A  careful  ••valuation  of 
the  low-lying  excited  states  of  U0i.  +  ,  particularly  those  charge  transfer 
states  that  arise  by  promotion  of  an  electron  from  (mainly)  ligand  MO's  to  a 
central  uranium  atom  MO,  indicates  that  such  transitions  lie  well  outside  of 
the  region  for  efficient  solar  pumping  of  U0^  +  .  However,  the  density  of 


electronic  states  for  U0(,+  lying  below  8  eV  is  immense,  and  additional  studies 


are  needed  for  this  system. 

A  preliminary  analysis  of  the  titanium/oxygen  system  was  also  carried  out 
with  emphasis  on  LWIR  for  TiO+.  An  unexpected  result  of  this  work  is  our 
conclusion  that  TiO+  is  a  weak  IR  radiator  with  fl0  =  4.4  x  10“ This  system 
is  thus  down  by  perhaps  a  factor  of  10  from  the  emission  intensity  of  the  U0+ 
and  UO  +  species.  However,  we  find  that  TiO<+  may  be  an  important  radiator 
and  further  character izat ion  of  this  system  is  required. 

Other  studies  within  DNA  chemistry  requirements  include  a  preliminary 
analysis  of  dielectronic  recombination  of  e  +  U+9  and  e  +  U0+ ,  which  are 
possible  charge  neutralization  processes,  and  an  analysis  of  the  kinetics  of 


the  N  +  0^  and  0  +  reactions  in  producing  vibrat  ional ly  hot  NO  molecules 
radiating  in  the  IR.  Our  analysis  of  dielectronic  recombination  indicates 
that  e  +  U+*  is  not  rate  competitive  with  possible  ion-neutral  charge  exchange 
reactions  but  that  e  +  U+  may  exhibit  a  large  kinetic  rate  for  neutral iza- 
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SECTION  1 


INTRODUCTION 

The  release  of  certain  chemical  species  into  the  upper  atmosphere  results 
in  luminous  clouds  that  display  the  resonance  e lec t r on ic-v i brat  ion-rot  at i on 
spectra  of  the  chemically  reacting  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies,  upon  re-entry  into  the 
atmosphere  of  artificial  satellites  and  missiles,  and  as  a  result  of  energv 
deposition  in  the  atmosphere  caused  by  nuclear  weapons  effects.  Of  particular 
interest  in  this  connection  is  the  observed  spectra  of  certain  metallic 
oxides.  Ffom  band  intensity  distributions  of  the  spectra,  and  knowledge  of 
the  f-values  for  electronic  and  vibrational  transitions,  the  local  conditions 
of  the  atmosphere  can  be  determined  (Reference  1).  Such  data  are  fundamental 
for  the  analysis  of  detection  and  discrimination  problems. 

Present  theoretical  efforts,  which  are  directed  toward  a  more  complete 
and  realistic  analysis  of  the  transport  equations  governing  atmospheric  relax¬ 
ation  and  the  propagation  of  artificial  disturbances,  require  detailed  infor¬ 
mation  of  thermal  opacities  and  LWIR  absorption  in  region  of  temperature  and 
pressure  where  both  atomic  and  molecular  effects  are  important  (References  2 
and  3).  Although  various  experimental  techniques  have  been  employed  for  both 
atomic  and  molecular  systems,  theoretical  studies  have  been  largely  confined 
to  an  analysis  of  the  properties  (bound-bound,  bound-free  and  free-free)  of 
atomic  systems  (References  4  and  3).  This  has  been  due  in  large  part  to  the 
unavailability  of  reliable  wavefunct ions  for  diatomic  molecular  systems,  and 
particularly  for  excited  states  or  states  of  open-shell  structure.  Only 
recently  (References  6-8)  have  reliable  procedures  been  prescribed  for  such 
systems  which  have  resulted  in  the  development  of  practical  computational 
programs . 

The  application  of  these  computational  methods  to  studies  of  the  elec¬ 
tronic  structure  and  radiation  characteristics  of  metal  oxides  has  been 
reported  for  several  of  the  lighter  systems  (References  9-11).  A  preliminary 
study  of  the  uranium/oxygen  system  has  been  reported  by  Michels  (Reference 
12)  which  identified  a  large  number  of  low-lying  molecular  states  for  both  the 
UO  and  (J0+  systems.  Of  particular  interest  was  the  discovery  of  two  struc¬ 
tures  for  U0+  that  resulted  from  two  different  spin-couplings  of  the  uranium 


•Si 


valence  electrons.  These  results  suggested  strong  LWIR  radiation  in  U0+  aris¬ 
ing  from  different  electronic  state  transitions. 

An  inherent  uncertainty  in  these  preliminary  calculations  was  present, 
owing  to  the  neglect  of  relativistic  effects  that  were  much  too  difficult  to 
include  in  molecular  calculations  at  that  period  of  time.  The  7s  valence 
electron  of  uranium,  and  its  corresponding  (J-bonding  molecular  orbital,  are 
highly  relativistic  in  nature  which  results  in  a  contracted  charge  density 
relative  to  that  which  would  occur  in  lighter  molecular  systems.  The  effect 
of  this  contraction  on  the  relative  positions  of  the  low-lying  electronic 
states  of  the  uranium/oxygen  system  can  now  be  calculated  with  some  degree  of 
confidence  using  newly  developed  relativistic  computer  codes. 

Because  of  inherent  difficulties  in  the  experimental  determination  of  the 
spectroscopy,  transition  probabilities  and  LWIR  radiation  for  metal  oxide 
systems  and  in  light  of  the  aforementioned  recent  progress  in  the  calculation 
of  relativistic  electronic  wave func t ions ,  especially  for  diatomic  systems,  a 
technical  program  for  calculating  these  properties  was  undertaken  for  the 
Defense  Nuclear  Agency  under  Contracts  DNAOO 1 -85 -C-0 1 20 ,  DNA00 1 -83 -00044  and 
DNAOO 1-82-C-001 5 .  The  emphasis  in  this  work  was  on  the  ions  of  uranium  and 
uranium  oxide  since  these  species  have  been  determined  to  be  important  radia¬ 
tors  in  the  LWIR  region.  These  studies  indicated  that,  in  addition  to  U0+  and 
U0^ + ,  the  doubly  ionized  species,  U++,  U0++  and  U0,  +  +,  should  also  be  consid¬ 
ered  because  of  their  role  in  charge  neutralization  processes  and  their  poten¬ 
tial  as  early-time  radiators. 

A  summary  of  the  results  obtained  to  date  for  the  uranium/oxygen  system 
indicates  that  both  U0+  and  UCh+  are  very  strong  LWIR  radiators  in  the  10-14 
Urn  region.  In  addition,  U()+ ,  which  is  a  potential  late  time  radiator  exhibits 
strong  visible  bands  and  should  be  efficiently  solar  pumped.  A  careful  evalu- 
at  ion  of  the  low-lying  excit’d  states  of  UO.+,  particularly  those  charge 
transfer  states  that  arise  by  promotion  of  an  electron  from  (mainly)  ligand 
MO's  to  a  central  uranium  atom  MO,  indicates  that  such  transitions  lie  well 
outside  of  the  region  for  etticient  solar  pumping.  However,  the  density  of 
electronic  states  tor  HO  +  lying  below  H  eV  is  immense,  and  additional  studies 
ire  needed  for  this  system. 

A  preliminary  analysis  ot  the  t i t an i um/oxvgen  system  was  also  carried  out 
with  emphasis  on  LWIR  tor  TiO+ .  An  unexpected  result  of  this  work  is  ,uir 


conclusion  that  TiO  is  a  weak  1R  radiator  with  f  t  ^  =  4.4  x  10  .  This  system 
is  thus  down  by  perhaps  a  factor  of  10  from  the  emission  intensity  of  the  U0+ 
and  UO/  species.  However,  we  find  that  TiOz+  may  be  an  important  radiator 
and  further  characterization  of  this  system  is  required. 

Other  studies  within  DNA  chemistry  requirements  include  an  analysis  of 
dielectronic  recombination  of  e  +  U+  and  e  +  U0+  which,  are  possible  charge 
neutralization  processes,  and  an  analysis  of  the  kinetics  of  the  N  +  0,  and 
0  +  N^,  reactions  in  producing  vibrat  ional  ly  hot  NO  molecules  radiating  in  the 
IR.  This  report  addresses  theoretical  studies  in  all  of  the  areas  described 
above . 

The  general  composition  of  this  report  is  as  follows.  In  Section  2,  we 
present  a  description  of  the  mathematical  methods  which  were  employed  in  this 
research.  Included  in  Section  2  are  sub-sections  which  deal  with  the 
construction  of  electronic  wavefunctions ,  the  calculations  of  expectation 
properties,  the  evaluation  of  molecular  transition  probabilities,  and  the 
calculation  of  electronic  wavefunctions  using  both  the  ab  initio  and  density 
functional  methods.  This  is  followed  by  Section  3  which  describes  the  inclu¬ 
sion  of  relativistic  effects  into  the  density  functional  (Xa)  method.  The 
calculated  results  and  pertinent  discussion  are  presented  in  Section  4. 
Recommendations  are  presented  in  Section  5. 


SECTION  2 


METHOD  OF  APPROACH  -  NONRF.LAT L V  1ST  1C  METHODS 


2.1  QUANTUM  MECHANICAL  CALCULATIONS. 


Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical 
calculations  which  must  be  carried  out  for  tin*  atomic  and  molecular  species. 

In  particular,  the  methodology  for  computing  the  minimum  energy  stable  struc¬ 
tures  must  be  carefully  analyzed.  For  added  clarity,  various  aspects  of  these 


calculations  are  discussed  below. 


Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy 
of  a  minimum  basis  set  for  constructing  quantitatively  correct  molecular  wave- 


functions  (References  13  and  14).  This  means  inner-shell  and  valence-shell 


orbitals  of  quantum  numbers  appropriate  to  the  atoms  (Is,  2s,  2p ,  for  C,  N,  0; 
etc.).  The  main  deficiency  of  the  minimum  basis  set  is  its  inability  to 
properly  describe  polarization  and  the  change  of  orbital  shape  for  systems 
which  exhibit  large  charge  transfer  effects.  Values  of  tlut  screening  param¬ 
eters  for  each  orbital  can  either  be  set  from  atomic  studies  or  optimized  in 
the  molecule;  the  latter  approach  is  indicated  for  studies  of  higher  preci¬ 
sion.  When  high  chemical  accuracy  is  required,  as  for  detailed  studies  of  the 
ground  or  a  particular  excited  state  of  a  system,  a  more  extended  basis  must 
be  used.  Double-zeta  plus  polarization  functions  or  optimized  MO's  are 
usually  required  as  a  minimum  representation  for  reliable  calculated  results 


of  chemical  accuracy. 


The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximal lv 
flexible  manner.  This  implies  the  construction  of  perturbation  expansions  or 
the  use  ot  Cl  wave  1  unc  t  i  oils  with  all  kinds  of  possible  orbital  occupancies,  so 
that  the  correlation  of  electrons  into  overall  states  can  adjust  to  an  opt  imum 


lorm  at  each  geometrical  conformation  and  for  each  state.  Except  when  we  1  1 - 
defined  pairings  exist,  as  for  closed  shell  and  exchange  dominated  systems,  a 


s t ng 1  e-con f  igu r at  i on  study  (even  ot  Hart ree-Fock  quality)  will  be  inadequate. 


Proper  electronic  states  tor  systems  composed  ot  1  i  gh  t  atoms  should 
possess  do  t  i  n  i  t  e  ••  i  , ’.onva  lues  o  1  the  spin  operal  or  S  as  well  as  an  appro 


ipr  i  at  •* 


m*  ■  t  r  i  c  a  I  svmmet  rv.  The  goomot  rical  svmmot  rv  can  he  coni  rol  led  hv  t  h* 


i  s  s  i  *  nine  nl  • )  !  orbitals  to  each  configuration,  but  the  spin  st  at  ••  rn  ist  h< 


obtained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (Reference  15)  for  projected  construction  of  spin  states  from 
orthogonal  orbitals,  and  programs  implementing  these  formulas  have  been  in 
routine  use  at  UTRC  for  several  years.  One  of  the  least  widely  appreciated 
aspects  of  the  spin-projection  problem  is  that  the  same  set  of  occupied 
spatial  orbitals  can  sometimes  be  coupled  to  give  more  than  one  overall  state 
of  a  given  spin  quantum  number.  It  is  necessary  to  include  in  calculations 
all  such  spin  couplings,  as  the  optimum  coupling  will  continuously  change  with 
changes  in  the  molecular  conformation.  This  is  especially  important  in 
describing  degenerate  or  near-degenerate  excited  electronic  states. 

In  the  sections  below,  we  describe  the  several  mathematical  approaches 
that  are  applicable  to  calculation  of  the  potential  energy  surfaces  for  chemi¬ 
cally  reacting  systems  and  to  the  subsequent  calculation  of  transition 
probabilities  and  f-numbers.  Since  several  different  approaches  are 
indicated,  we  describe  their  expected  regions  of  applicability. 


2.1.1  Method  of  Ab  Init io  Calculation. 

2. 1.1.1  Born-Oppenheimer  Separation.  For  a  system  of  n  electrons  and  N 
nuclei,  and  considering  only  electro-static  interactions  between  the  parti¬ 
cles,  we  have  for  the  total  Hamiltonian 


where 
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and  where  m  ,  m  ,  M  ,  are  the  masses  of  the  electron,  atom  «  and  combined 
e  a  T 

system  mass,  respectively.  Now  since  the  ratios  m^/m^  and  m^/M^,  are  both 
small,  (2  •  L0_o  -  5  •  10~4)  we  can  effect  a  separation  of  the  electronic  and 
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JV-v; 


nuclear  coordinates  treating  the  total  wavefunction  as  a  product  of  a  nuclear 
and  an  electronic  part.  We  have 


k 


(3) 


where  b  (r  ,  R  )  is  an  electronic  wavefunction  parametric  in  the  nuclear 
k  ~n  ~N 

coordinates  as  given  in  Equation  (3)  and  X,  (R„)  are  nuclear  motion  wavefunc- 

o  t  k  ~N 

tions  which  satisfy  (neglecting  terms  of  the  order  of  me/Ma) 
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The  cross  term  in  can  be  eliminated  by  a  proper  change  of  vari¬ 

ables  and  Equation  (4)  then  reduces  to  a  3N-3  dimensional  Schrodinger 
equat ion. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow 

relative  to  the  electron  velocity,  this  decoupling  of  electron  and  nuclear 

motion  is  valid  and  is  referred  to  as  the  adiabatic  approximation.  Equation 

(3)  thus  defines  an  electronic  eigenstate  b. (r  ,  R  ),  parametric  in  the 

k  ~n  ~N 

nuclear  coordinates,  and  a  corresponding  eigenvalue  E,  (R  )  which  is  taken  to 

k  ~N 

represent  the  potential  energy  curve  or  surface  corresponding  to  state  k. 


2 . 1 .  1 . 2  Variational  Methods.  By  an  jit>  initio  method  is  meant  one  that 
starts  from  a  zero-order  Hamiltonian  which  is  exact  except  for  relativistic 
and  magnetic  effects,  and  which  involves  the  evaluation  of  electronic  energies 
and  other  relevant  quantities  for  wave funct ions  which  are  properly 
antisymmetrized  in  the  coordinates  of  all  the  electrons.  For  a  system 
containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian  depends 
parametrically  on  the  nuclear  positions  and  is  of  the  form 


^  y  ,2_£  £  J  +  £  1  _1  +  £  ..  1 

2  i  =  i  i=i  j  =  i  r~,  -  fTj  i<i<j  R,  -R,  '<!<)  |  r ,  -  r j 


(5) 


where  and  R^  are  the  charge  and  position  of  nucleus  i,  r.  is  the  position 
of  electron  j,  and  Vj^  is  the  Laplacian  operator  for  electron  j.  All 
quantities  are  in  atomic  units,  i.e.  lengths  in  bohrs ,  energies  in  hartrees 
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(1  hartree  -  2  Rydbergs).  The  many-electron  wave func t ion  consists  of  one,  or 
a  linear  combination,  ^  =  i  c  ,  of  terms  of  the  form 

II  ”  " 


where  each  <f> ^  .  is  a  spatial  orbital,  $  is  the  antisyrametr izing  operator,  0  is 
the  spin-projection  operator  for  spin  quantum  number  S,  and  9  is  a  product 
of  a  and  8  one-electron  spin  functions  of  magnetic  quantum  number  Ms .  No 
requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital,  so 
linear  combinations  of  the  form  given  by  Equation  (6)  can  describe  a 
completely  general  wavefunction.  The  spatial  orbitals  $  ^  may  be  whatever 
basis  orbitals  have  been  introduced,  arbitrary  linear  combinations  thereof,  or 
specific  linear  combinations  determined  pursuant  to  the  particular  calcula- 
tional  method  in  use. 

The  spatial  orbitals  <p  .  ,  the  spin  functions  9  and  the  coefficients 

U  l  pm  ’ 

s 

of  different  'J'p,  if  a  linear  combination  of  is  used,  may  be  explicitly 
determined  by  invoking  the  variational  principle.  Various  specific  methods 
are  described  below  for  determining  wave func t ions .  However,  we  should  first 
observe  that  the  adequacy  of  an  ab  initio  calculation,  or  for  that  matter  any 
energy  calculation,  will  depend  crucially  upon  the  extent  to  which  the  wave- 
function  can  be  qualitatively  appropriate.  Some  of  the  considerations 
surrounding  the  choice  of  a  wavefunction  are  the  following: 

(i)  necessity  that  the  wavefunction  possess  sufficient  flexibility  to  be 
able  to  describe  dissociation  to  the  correct  atomic  and  molecular  fragments  as 
various  internuclear  separations  are  increased; 

(ii)  maintenance  of  equivalent  quality  of  calculation  for  nuclear  geome¬ 
tries  differing  in  the  nature  or  number  of  chemical  bonds; 

(iii)  ability  to  describe  degenerate  or  near-degenerate  electronic  states 
when  they  are  pertinent; 

(iv)  ability  to  describe  different  electronic  states  to  equivalent 
accuracy  when  their  interrelation  (e.g.,  crossing)  is  relevant,  in  particular, 
ability  to  describe  ionic-valence  state  mixing; 

(v)  ability  to  represent  changes  in  the  coupling  of  electron  spins  as 
bonds  are  broken  or  reformed. 


The  foregoing  considerations  indicate  that  it  will  often  be  necessary  to 
consider  wavefunctions  with  more  than  a  minimum  number  of  singly-occupied 
spatial  orbitals,  and  that  there  will  be  many  potential  curves  or  surfaces  for 
which  a  wavefunction  consisting  of  a  single  cannot  suffice.  It  will  then 
be  necessary  to  allow  mixing  of  'f'p  with  different  degrees  of  orbital  spatial 
occupancy  so  as  to  obtain  smooth  transitions  from  the  occupancies  characteris¬ 
tic  of  separated  atoms  or  molecules  (or  ions)  to  those  characteristic  of  a 
compound  system  or  a  different  fragmentation. 

Another  implication  of  the  considerations  surrounding  the  choice  of  a 
wavefunction  is  related  to  the  treatment  of  electron  spin.  Not  only  is  it 

*  •  9 

necessary  to  require  that  the  wavefunction  be  an  eigenfunction  of  S  and  Sz 
but  it  is  also  necessary  to  take  account  of  the  fact  that  under  many  condi¬ 
tions,  there  will  be  more  than  one  spin  eigenfunction  of  given  S  and  Ms . 

The  different  spin  eigenfunctions  correspond  to  different  couplings  among  the 
individual  spins.  Since  reactive  processes  involve  the  breaking  and  forming 
of  electron-pair  bonds,  they  must  necessarily  be  accompanied  by  reorganiza¬ 
tions  of  the  spin  coupling.  A  failure  to  take  account  of  this  will  lead  to 
qualitatively  inappropriate  wavefunctions. 

In  Hartree-Fock  calculations  ^(R)  is  restricted  to  a  single  which  is 

u 

assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 

orbitals  <J>  .  are  then  selected  to  be  the  linear  combinations  of  basis  orbitals 

u  1 

best  satisfying  the  variational  principle.  Writing  ^  •  =  E  a  .  X  ,  the  a  . 

J  0  r  r  &  u  l  V  l  V  Vl 

are  determined  by  solving  the  matrix  Hartree-Fock  equations 


?  sx„a„,  (eac*  M 


where  is  the  orbital  energy  of  ^ . 

The  Fock  operator  F^y  has  been  thoroughly  discussed  in  the  literature  (Refer¬ 
ence  16)  and  depends  upon  one-  and  two-electron  molecular  integrals  and  upon 
the  a  ..  This  makes  Equation  (7)  nonlinear  and  it  is  therefore  solved  itera- 

*  L 

tively.  UTRC  has  developed  programs  for  solving  Equation  (7)  for  both  closed 
and  open-shell  systems,  using  basis  sets  consisting  of  Slater-type  atomic 
orbitals.  Examples  of  their  use  are  in  the  literature  (Reference  7). 


AV-V-AA/-.- 
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In  configuration  interaction  calculat ions ,  the  overall  wavefunction  has 

than  one  term,  0  ,  and  the  c  are  de 

U  u 

principle  to  obtain  the  secular  equation 


more  than  one  term,  0^,  and  the  c^  are  determined  by  invoking  the  variational 


X  ( Hfj.v  ~  )  c„  =  0  (each /x) 


(8) 


where 


=  / VjiMJ  (R)^v(R)^r 

S J (R)Vv  (R)  dr 


(9) 


Equation  (8)  is  solved  by  matrix  diagonalization  using  either  a  modified 
Givens  method  (Reference  17)  or  a  method  due  to  Shavitt  (Reference  18)  or 
Raffenetti  (Reference  19). 

The  matrix  elements  H  and  may  be  reduced  by  appropriate  operator 

algebra  to  the  forms 
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where  P  is  a  permutation  and  e  its  parity.  The  sum  is  over  all  permutations. 

P 

<0  1(9  P|9  >  is  a  "Sanibel  coefficient"  and  the  remaining  factors  are  spatial 

M1  s  M 

integrals  which  can  be  factored  into  one-  and  two-electron  integrals.  If  the 

<t>  .are  orthonormal,  Equations  (10)  and  (11)  become  more  tractable  and  the  H 

111  u  v 

and  may  be  evaluated  by  explicit  methods  given  in  the  literature 

(Reference  15).  Computer  programs  have  been  developed  for  carrying  out  this 
procedure,  and  they  have  been  used  for  problems  containing  up  to  106  total 
electrons,  10  unpaired  electrons,  and  several  thousand  configurations. 


The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set 

of  <(>  .  for  which  the  molecular  integrals  can  be  calculated.  Programs  devel- 

M  l 

oped  by  UTRC  make  specific  provision  for  the  choice  of  the  as  Slater-type 

atomic  orbitals,  as  Gaussian-type  orbitals,  as  symmetry  molecular  orbitals,  as 
Hartree-Fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method 
of  calculation  are  evaluated  for  STO's  by  methods  developed  by  this  Center 
(Reference  20).  For  Gaussian  orbitals,  either  the  Carnegie-Me l Ion  integral 
package  (Reference  21)  or  the  integral  routines  incorporated  in  the  GAMES S 
program  (Reference  22)  can  be  employed.  All  needed  computer  programs  have 
been  fully  tested  at  UTRC. 


2. 1.1.3  Configuration  Selection.  Using  a  double-zeta  plus  polarization 
basis  set  of  one-electron  functions,  a  typical  system  can  have  of  the  order  of 
10J  configurations  in  full  Cl  (that  resulting  from  all  possible  orbital 
occupancies).  It  is  therefore  essential  to  identify  and  use  the 
configurations  describing  the  significant  part  of  the  wave func t ion .  There  are 
several  ways  to  accomplish  this  objective.  First,  one  may  screen  atomic- 
orbital  occupancies  to  eliminate  conf igurat ions  with  excessive  formal  charge. 
Alternatively,  in  a  molecular-orbital  framework,  one  may  eliminate 
configurations  with  excessive  numbers  of  anti-bonding  orbitals.  A  third 
possibility  is  to  carry  out  an  initial  screening  of  configurations,  rejecting 
those  whose  diagonal  energies  and  interaction  matrix  elements  do  not  satisfy 
energy  significance  criteria. 

Another  common  method  of  classifying  configurations  is  to  examine  the 
total  number  of  orbitals  in  the  wavefunction  that  differ  from  the  SCF 
reference  wavefunction.  We  write  the  Cl  wavefunction  as 
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where  a  single  excitation  function,  <f>  .  ,  differs  from  the  SCF  reference  <!'  by 

1  i)  0 

one  orbital  and  the  double  excitation  function  •!>  .  by  two,  and  so  on.  Hie  Cl 

.1 

coefficients,  C{,  are  then  determined  var  iat  ional ly  to  yield  the  lowest 
possible  total  energy.  In  the  limit  of  a  complete  basis  set  ( N  *-"° )  ,  and  where 
all  possible  substitutions  are  included  in  <1^,  ,  the  variational  energy 

approaches  the  correct  non  re  1  a t  i  v i s t i c  Horn -Oppenhe ime r  result. 


Errors  arise 


from  a  truncation  of  the  functions  used  to  determine  the  SCF  reference  wave- 
function  and  from  the  truncation  of  the  excitation  functions  series  in  1J1 

C I 

The  reference  wavefunction  1J1  will  typically  be  the  same  for  a  Cl  or  Many  Body 

o 

Perturbation  Theory  (MBPT)  calculation;  however,  for  a  molecule  even  as  small 
as  water,  ip  becomes  a  function  of  a  very  large  number  of  basis  functions. 

L*  1 

Because  of  this,  one  must  truncate  the  expansion  and  eliminate  unimportant 
configurations.  In  general  one  can  show  that 


<*oWI*T>  =  <K^Q>  =  0 


In  words,  the  matrix  elements  between  the  reference  wavefunction  and  triple 
and  quadruple  excitations  is  zero.  Thus  to  first  order,  only  single  and 
double  excitations  contribute.  Because  of  this,  many  Cl  calculations  attempt 
to  include  all  single  and  double  excitations  in  the  expressions  for  ij/  .  To 
go  beyond  this,  generally  more  than  one  reference  wavefunction  is  used. 
Programs  to  handle  configurations  on  all  the  above  criteria  are  available  at 
UTRC . 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Reference  23)  or  multiconfigura¬ 
tion  SCF  (Reference  24)  concepts.  To  implement  the  natural-orbital  approach, 
an  initial  limited-CI  wavefunction  is  transformed  to  natural-orbital  form,  and 
the  resulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  desired  result 
is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller  number  of 
significant  terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome, 
but  in  principle  more  effective.  It  yields  the  optimum  orbital  choice  for  a 
preselected  set  of  configurations.  This  approach  works  well  when  a  small 
number  of  dominant  configurations  can  be  readily  identified.  The  method  is 
described  briefly  below. 

2. 1.1.4  Multiconfieuration  -  Self  Consistent  Field  Method  (MC-SCF).  The 


Hartree-Fock  self  consistent  field  method  has  been  proven  to  be  a  powerful 
tool  for  the  calculation  and  understanding  of  many  ground  state  properties  of 
molecules  in  the  vicinity  of  their  equilibrium  structure.  However,  in  most 
cases  the  one  determinant  Hartree-Fock  approximation  is  not  adequate  to 
properly  describe  the  dissociation  of  molecular  bonds.  Also,  many  excited 
states  cannot  be  represented  by  a  single  configuration  wave  function.  In 
order  to  calculate  properties  of  such  states  or  to  investigate  the  formation 
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of  molecular  bonds  one  often  needs  multiconfiguration  wave  functions  for  which 
both  the  linear  coefficients,  of  the  configuration  expansion 
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as  well  as  the  set  of  orthonormal  molecular  orbitals  { d> j  } »  from  which  the 
configurations  are  constructed,  are  optimized  according  to  the  varia¬ 

tional  principle.  As  is  well  known,  this  "MC-SCF"  problem  presents  many  more 
difficulties  than  the  simple  one  determinant  Hartree-Fock  case  and  much  work 
has  been  devoted  to  obtaining  convergent  solutions  during  the  last  decade. 

The  difficulties  mainly  arise  from  the  fact  that  for  general  MC-SCF  wave 
functions  the  energy  is  not  invariant  with  respect  to  rotations  between 
occupied  orbitals.  Hence,  instead  of  a  relatively  simple  pseudo-eigenvalue 
equation  in  the  one  determinant  case,  the  set  of  coupled  Fock  equations 
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with  the  hermiticity  conditions 
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has  to  be  solved.  £ . .  are  Lagrange  multipliers  which  account  for  the  ortho¬ 
normality  constraints  imposed  on  the  orbitals.  The  Fock  operators  F„  depend 
on  the  orbitals  { 4>  £ }  and  the  set  of  Cl  coefficients  { C  £  } .  In  analogy  to 
the  one  determinant  case  many  attempts  have  been  made  to  solve  these  equations 
iteratively  by  keeping  the  Fock  operators  fixed  in  each  iteraction  step.  Then 
the  Lagrange  multipliers  can  be  expressed  by  coupling  operators  constructed 
such  that  the  Fock  equations  are  transformed  into  pseudo-eigenvalue  equations 
yielding  the  improved  orbitals.  These  are  used  in  a  second  step  to  determine 
new  CL  coefficients  by  diagonalizing  the  Cl  matrix  <<l>  ^  >  •  The  convergence 

of  these  algorithms,  however,  has  often  been  found  to  be  poor. 

A  cprmnd  group  of  MC-SCF  methods  is  based  on  the  generalized  Brillou’n 
theorem.  In  these  methods,  the  orbital  changes  are  derived  from  the  coeffi¬ 
cients  of  a  Cl  expansion  consisting  of  the  MC-SCF  wavefunction  and  all  one- 
electron  singly  excited  configurations.  A  computer  program  (ALLS)  imple¬ 
menting  this  method  has  been  developed  by  Ruedenberg,  et  al  (Reference  25).  A 
somewhat  more  elegant  program  (GAMESS)  which  utilizes  analytical  gradients  is 
also  available  (Reference  22).  More  recently,  various  methods  have  been 
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proposed  which  are  based  on  direct  minimization  of  the  energy,  avoiding  the 
Fock  operators  altogether.  Such  methods  are  now  being  studied  in  our  labora¬ 
tory  and  will  be  incorporated  into  our  existing  computer  programs  if  they 
prove  to  be  highly  efficient. 

2. 1.1. 5  Many-Body  Perturbation  Theory  (MBPT).  In  MBPT  we  again  begin  with 
the  SCF  wavefunction  (4^)  as  our  reference  and  attempt  to  account  for  Ec , 
the  correlation  energy.  The  concept  of  excitation  functions  described  in  the 
above  section  on  Cl  calculations  carries  over  to  MBPT  calculations.  In  MBPT 
one  can  write  the  wavefunction,  4>p,  as 

«  eT |  <j>Q>  (17) 

where  T  is  an  excitation  operator  defined  as 

T  =  Ts  +  td  +  TT  +  TQ  •  •  •  (18) 

wnere  S,  D,  T,  and  Q  refer  to  single,  double,  triple,  and  quadruple  substi¬ 
tutions  respectively.  One  can  write  Tn ,  in  general,  where  n  refers  to  the 
number  of  excitations  as 


T  =  — 
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where  a,b,c,...  are  excited  orbitals  and  i,j,k,...  are  orbitals  occupied  in 
<f>0.  The  total  energy  is  now  given  by 


EMBPT  <<J)o  eT  I  <j5o> 


To  evaluate  Ewnr>_,  the  t.’’  from  above  must  be  determined.  An  equivalent 
Mori  1  .  .  . 

expression  that  makes  the  perturbation  expansion  clearer  is 

00 

EMBPT  =  '•  ~^o)  ^]k|4>0> 


where  the  sum  is  over  only  so  called  linked  diagrams,^  has  eigenfunctions 
4>  ,  and  the  expansion  is  of  orders  in  the  perturbation,  V  =  (cM -  G).  The 
k  =  0  term  gives  the  reference  energy  and  for  k  >  0  correlation  corrections 


IJF*. 


are  included.  In  practice  the  MBPT  totaL  energy  is  calculated  by  truncating 

T 

the  T  operator  expansion  and  projecting  de  n|<}>^>  onto  the  appropriate  n- 
space.  This  leads  to  a  set  of  nonlinear  coupled  equations  for  the  t3’’’ 
coefficients  which  correspond  to  the  Cl  expansion  coefficients.  The  equations 
are  solved  iteratively  and  evaluated.  In  practice  T4  is  an  upper  limit 

that  corresponds  to  quadruple  substitutions.  The  series  is  an  oscillatory 
convergent  sum,  which  in  practice  has  proven  to  be  at  least  as  accurate  as 
E(j[  with  single  and  double  excitations  included. 

The  best  possible  method  to  use  would  be  a  full  Cl  (all  possible  configu¬ 
rations)  with  a  complete  basis  set.  However,  since  the  number  of  configura¬ 
tions  is  proportional  to  <n  |  ^  ,  where  n  is  the  number  of  basis  functions  and 
£  is  the  level  of  excitation,  it  is  usually  prohibitive  to  even  include  all 
single  and  double  excitations.  This  truncation  causes  the  loss  of  size 
consistency,  which  implies  that  the  energy  calculated  for  A  and  B  as  a 
molecular  system,  but  far  apart,  is  the  sum  of  the  energy  calculated  for  A  and 
B  separately.  In  a  size  extensive  calculation  the  energy  is  proportional  to 
the  size  of  the  system.  These  properties  are  very  important  if  one  wishes  to 
compute  correct  relative  energies  on  a  potential  energy  surface,  a  necessary 
criteria  for  defining  the  reactive  charge  transfer  pathways  of  interest  in 
this  research  program.  MBPT,  on  the  other  hand,  is  guaranteed  to  have  the 
correct  size-dependence  because  the  expansions  contain  only  the  so-called 
linked  diagrams.  In  addition,  because  of  the  computational  efficiency  of 
MBPT,  calculations  can  be  performed  up  to  fourth  order  in  the  perturbation 
expansion  and  include  single,  double,  triple,  and  quadruple  excitations  in  the 
calculation  of  the  correlation  energy.  Such  calculations  are  often  performed 
with  a  split  valence  plus  polarization  basis  set.  For  the  light  element 
compounds  which  exhibit  ionic  bonding,  the  inclusion  of  diffuse  basis  func¬ 
tions  and  possibly  higher  polarization  d  or  f-functions  may  be  required  for  a 
quantitative  treatment.  Again,  all  of  the  necessary  programs  to  handle 
expansions  up  to  MPA  (SI)TQ),  with  inclusion  of  up  to  f-functions  are  opera¬ 
tional  and  in  routine  use  at  UTKC. 

2.1.2  Density  Functional  Approach. 

The  X«  method  (Reference  26)  for  the  electronic  structure  of  atoms, 
molecules,  clusters  and  solids  is  a  local  potential  model  obtained  by  making  a 
simple  approximation  to  the  exchange  -  correlation  energy.  Although  this 
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method,  as  currently  implemented,  is  less  accurate  than  standard  at>  initio 
approaches,  it  is  very  useful  for  heavy-atom  molecular  systems.  For  such 
systems  (Z  >  30),  relativistic  effects  become  increasingly  important  and  their 
incorporation  into  density  functional  methods  is  much  easier  than  into  ab 
init io  methods.  If  we  assume  a  nonre 1  at iv is t  ic  Hamiltonian  with  only  electro¬ 
static  interactions,  it  can  be  shown  that  the  total  energy  E  of  a  system  can 
be  written  exactly  (Reference  24)  (in  atomic  units)  as 


I  I  fJ.  I /A  I  1  /x^l/  r 


*  T- 1  ni"i<u,  hi  tJu‘ui>  +  E’ 


This  expression  is  exact  provided  the  u^  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  Equation  (22)  represents  the  kinetic 
and  electron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion 
energy.  The  sums  (u,v)  are  over  all  the  nuclear  charges  in  the  system.  The 
third  term  is  the  electron-electron  repulsion  term,  which  represents  the  clas¬ 
sical  electrostatic  energy  of  the  charge  density  0  interacting  with  itself, 


where 


p(i)  =Z  n.u*(D  UjO) 


The  last  term  E  represents  the  exchange  correlation  energy  and  can  be 


expressed  formally  as 


E.c  P«  ":?i  dr? 


where  p  (1,  2)  represents  the  exchange-correl at  ion  hole  around  an  electron  at 
xc 

position  1.  In  the  exact  expression,  p  is  dependent  on  the  second-order 

density  matrix.  In  the  Hartree-Fock.  approximation  E  is  the  exchange  energy, 

P  represents  the  Fermi  hole  due  to  the  exclusion  principle  and  depends  only 
xc 


Em 

$1 

r» 


i 


on  the  first-order  density  matrix.  [n  the  Xm  method,  we  make  a  simpler 

assumption  about  P  .  If  we  assume  that  the  exchange-corre Lat ion  hole  is 
xc 

centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown  that 
the  exchange-correlation  potential 


U  =  f ^ 

*c  J  ri2 


dr 


is  inversely  proportional  to  the  range  of  the  hole,  rs,  where  rs  is 
tie  fined  by 


3  r  s  P 


e3  />(')=! 


Therefore,  in  the  X(j  model,  the  potential  II  is  proportional  to  P  (r).  We 
define  a  scaling  parameter  a  such  that 


Uxa,^=  “  2  ( SpiD/eTr)1 


The  expression  in  Equation  (27)  is  defined  so  that  «  =  2/3  for  the  case  o f  a 
free  electron  gas  in  the  Hartree-Fock  model  (Reference  28)  and  <x  =  1  for  the 
potential  originally  suggested  by  Slater  (Reference  29).  A  convenient  way  to 
choose  this  parameter  for  molecular  and  solid  state  applications  is  to  opti¬ 
mize  the  solutions  to  the  X<j  equations  in  the  atomic  limit.  Schwarz 
(Reference  30)  has  done  this  for  atoms  from  z  =  l  to  z  =  41  and  found  values 
between  2/3  and  1. 

In  the  "spin  polarized"  version  of  the  X<t  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  onlv  with  a 
potential  determined  by  the  charge  density  of  the  same  s p i n .  In  this  case*  the 
contribut  ion  to  the  total  energy  is  summed  over  the  two  spins,  s  =  *  1/2. 


E.c  =  -r?/^i"u.a.s(""h 


V  .>  N 

VvVV 


where  the  potential  is  spin-dependent 


IL  (>)=- 

s 


9a 


(3/>(l)/4ir) 


1/3 


(29) 


and  Ps  is  the  charge  density  corresponding  to  electrons  of  spin  s.  The 
spin  polarized  Xa  model  is  useful  for  describing  atoms  and  molecules  with 
open-shell  configurations  and  crystals  which  are  ferromagnetic  or  anti¬ 
ferromagnetic  . 

Once  one  has  made  the  Xa  approximation  to  the  total  energy  functional  E 
in  Equation  (22),  then  the  rest  of  the  theory  follows  from  the  application  of 
the  variational  principle.  The  orbitals  ut  are  determined  by  demanding  that 
E  be  stationary  with  respect  to  variations  in  u{ .  This  leads  to  the  set  of 
one-electron  Xa  equations 


♦  I 

/* 


/ 


^rt* 
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uxaJ  ui  = 


u- 


(30) 


where  is  the  one-electron  eigenvalue  associated  with  u^.  Since  3(r)  is 
defined  in  terms  of  the  orbitals  u^,  Equation  (30)  must  be  solved  itera¬ 
tively,  until  self-consistency  is  achieved.  Empirically,  if  one  takes  as  an 
initial  guess  that  p  is  approximately  a  sum  of  perimposed  atomic  charge 
densities,  then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor 
of  2/3  multiplying  the  potential  is  a  result  of  the  linear  dependence  of  E 

xc 

on  P.  This  also  has  a  consequence  that  the  Xa  eigenvalues  do  not 
satisfy  Koopman's  theorem,  i.e.,  they  cannot  be  interpreted  as  ionization 
energies.  However,  it  can  be  shown  that  the  are  partial  derivatives  of 
the  total  expression  of  Equation  (22)  with  respect  to  the  occupation  number, 


t 


d  E 
dn. 


(31) 


If  E  were  a  linear  function  of  n^,  then  Koopmans'  theorem  would  hold. 
However,  because  of  the  dominant  Coulomb  term,  E  is  better  approximated  by  a 
quadratic  function  in  n{.  This  leads  to  the  "transition  state"  approxima¬ 
tion  which  allows  one  to  equate  the  difference  in  total  energy  between  the 
state  (n.,  n.)  and  (n.-  1,  n.+  1)  to  the  difference  in  the  one-electron 


energies  e.  -  £ .  calculated  in  the  state  (n.-  1/2,  n.+  1/2).  The  error  in 
J  t  L  J 

this  approximation  is  proportional  to  third-order  derivatives  of  E  with 

respect  to  n£  and  n j ,  which  are  usually  small  (Reference  31).  The  main 

advantage  of  using  the  transition  state  rather  than  directly  comparing  the 

total  energy  values  is  computational  convenience,  especially  if  the  total 

energies  are  large  numbers  and  the  difference  is  small. 

The  relationship  of  Eq .  (31)  also  implies  the  existence  of  a  "Fermi 
level"  for  the  ground  state.  This  can  be  seen  by  varying  E  with  respect  to  n. 

V  .  1 

under  the  condition  that  the  sum  £  n^  is  a  constant,  i.e., 

»[e-XI  »J  =  °  (32) 

implies  8E/3n£  =  X;  where  ^  is  a  Lagrangian  multiplier.  This  implies  that 
the  total  energy  is  stationary  when  all  the  one-electron  energies  are  equal. 
However,  the  occupation  numbers  are  also  subject  to  the  restriction 
0  <  °i  <  1.  This  leads  to  the  following  conditions  on  the  ground  state 
occupation  numbers: 

Ei  <  X*)*  n  ^  =  1 

e{  >  A*)*  nL  =  0  (33) 

e  ^  =  X  • )  •  0<n  £_<  1 

In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with  X 
representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n£  are  either  0  or  1 ,  the  Xq  model  pre¬ 
dicts,  in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level.  In 
particular,  this  will  occur  in  a  system  (such  as  transition  metal  or  actinide 
atom)  which  has  more  than  one  open  shell. 

The  Xq  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree-Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  Equation  (30)  is  orbital-independent  and  local,  i.e.,  it  is  the 
same  for  all  electrons  (except  in  the  spin-polarized  Xn  theory)  and  is  a 
multiplicative  operator.  On  the  other  hand,  the  Hartree-Fock  potential  is 


nonlocal,  or  equivalently,  there  is  a  different  local  potential  for  each 
orbital.  This  involves  a  great  deal  more  computational  effort,  especially  for 
systems  described  by  a  large  number  of  orbitals.  It  has  been  shown  (Reference 
32)  that  the  X&  orbitals  for  the  first  and  second  row  atoms  are  at  least  as 
accurate  as  a  double-zeta  basis  set,  and  are  probably  better  for  larger  atoms 
which  involve  electrons  with  £  _>  2. 

2.  The  orbital-independent  Xa  potential  leads  to  a  better  one-electron 
description  of  electronic  excitations  of  a  system.  Both  the  unoccupied 

(n .  =  0)  and  occupied  (n .  =  1)  eigenfunctions  are  under  the  influence  of  the 
same  potential  resulting  from  the  other  N-l  electrons.  The  Hartree-Fock  vir¬ 
tual  orbitals  see  a  potential  characteristic  of  the  N  occupied  orbitals,  and 
therefore  are  not  as  suitable  for  describing  the  excited  states.  Actually, 
although  the  ground  state  virtual  eigenvalues  are  usually  a  good  description 
of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
potential  where  one-half  an  electron  has  been  removed  from  the  system  gives  a 
much  better  first-order  picture  of  these  levels  (Reference  33). 

3.  As  has  been  shown  by  Slater  (Reference  34),  the  Xq  model  rigorously 
satisfies  both  the  virial  and  Hel lman-Feynman  theorems,  independent  of  the 
value  of  the  parameter  a.  This  is  convenient  for  calculating  the  force  on  a 
nucleus  directly  in  terms  of  a  three-dimensional  integral,  rather  than  the 
six-dimensional  integrals  in  the  expression  for  the  total  energy  of  Equation 
(22)  . 


2. 1.2.1  Computational  Aspects  of  the  X^  Method.  In  application  of  the  Xq 

model  to  finite  molecular  systems,  there  are  two  practical  aspects  of  the 
calculations  which  must  be  considered.  The  first  concerns  the  choice  of  the 
integration  framework  for  describing  the  molecular  wave funct ions  and  the 
second  deals  with  the  choice  of  the  exchange  parameter,  a,  in  different 
regions  of  space. 

In  computations  with  heteronuc lear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
spheres  of  integration  at  a  given  internuclear  separation,  the  degree  of 
sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic  spheres 
and  the  intersphere  region. 
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It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a  heteronuc Lear  diatomic  molecule  introduces  changes  in  the  total 
energy  that  can  be  large  on  a  chemical  scale  (~  1  eV) .  A  choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a  good 
estimate  for  these  calculations.  The  value  of  the  exchange  parameter,  o i,  and 
the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  calculations 
for  crystals  where  the  geometry  is  fixed.  However,  to  develop  a  potential 
curve,  the  molecule  description  needs  to  change  substantially  as  the  inter- 
nuclear  separation  varies  and  the  changing  sphere  radii  include  varying  frac¬ 
tions  of  the  total  molecular  charge  (Reference  35).  Studies  made  at  UTRC  have 
shown  that  at  any  given  separation  the  total  energy  calculated  from  the  Xa 
model  is  a  minimum  at  the  radii  ratio  where  the  spherically  averaged  poten¬ 
tials  from  the  two  atomic  centers  is  equal  at  the  sphere  radius. 

VrSl>  *  Vr,2>  <34) 

This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic 
species,  the  long  range  tail  of  the  potential  must  go  like  +2/R  for  one  ion 
and  -2/R  (in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separa¬ 
tions,  the  tails  of  the  potential  cannot  match  well.  However,  at  reasonable 
separations,  the  1 / R  character  of  the  potential  does  not  invalidate  the  poten¬ 
tial  match  criterion  for  radii  selection.  This  match  for  the  atomic  poten¬ 
tials  is  applied  to  the  self-consistent  potentials. 

In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an 
overlap  region  appears  in  the  vicinity  of  the  bond  (Reference  36).  Studies 
made  at  UTRC  show  that  the  contribution  to  the  total  molecular  energy  from  the 
exchange  integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere  over¬ 
lap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 
neutral  and  for  ionic  molecular  constituents.  However,  for  light  atoms,  the 
value  of  a  which  best  reproduces  Hartree-Fock  results  varies  subs t ant i a  1 1 v 
with  ionic  ity.  In  argon,  the  following  table  compares,  for  the  neutral  atom 


v  vi  v  vw  v.  ^^.vtn?oy. 


and  the  positive  ion,  the  HF  energy  and  the  X a  energy  calculated  for  several 
values  of  a. 


a 

Xa  Energy 

HF  Energy 

Ar  0 

.72177 

526.8176 

526.8173 

Ar+1/2 

.72177 

526.5857 

_ 

.72213 

526.6007 

- 

Ar+  1 

.72177 

526.2447 

_ 

.72213 

526.2596 

- 

.72249 

526.2745 

526.2743 

The  optimum 

value  of  a 

changes  even 

more  rapidly  in  the 

fluorine  atom  ; 

from  0.73732 

for  F  0  to 

.72991  for  F‘ 

').  Since  the  total 

energy  depends 

ly  on  a,  this  parameter  must  be  chosen  carefully. 

The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average 
of  the  atomic  exchange  parameters  from  the  two  constituents.  At  small  inter- 
nuclear  separations,  the  optimum  radius  for  an  atomic  sphere  frequently  places 
significant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that  is 
still  strongly  associated  with  its  original  center  rather  than  being  trans¬ 
ferred  to  the  other  center  or  associated  with  the  molecular  binding  region. 

To  best  account  for  these  cases  the  weighting  coefficients  are  chosen  to 
reflect  the  origin  of  the  charge  in  the  intersphere  (or  outersphere  region), 


<QSl  -  O  +  as.  (Qc  -  Q ,°> 


inter sphere 


(Q„  -  Q  )  +  (Q 


(35) 


where  (Q  -  Q° )  is  the  charge  lost  from  sphere  i  relative  to  its  atomic  value 
i 

(or  ionic  value)  Q.  and  a  is  the  atomic  exchange  parameter  for  sphere  i. 

i 

This  value  for  a.  ,  is  calculated  dynamically  -  it  is  updated  after 

mtersphere 

each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  wouLd  be 
small,  the  o's  for  small  atoms  vary  rapidly  with  7  (and  with  ionicitv).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
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calculated  for  the  molecule  but  also  in  some  cases  affects  the  distribution  of 
charge  between  the  atomic  spheres  and  the  intersphere  region. 


2.1.3  Ab  Initio  Gaussian  Wavefunct ions . 

Owing  to  the  complexity  of  evaluating  multicenter  electron  repulsion 
integrals  over  Slater-type  (exponential)  orbitals,  various  groups  have  adopted 
a  computational  approach  to  electronic  structure  calculations  based  on 
Gaussian  orbitals.  A  highly  developed  program,  named  Gaussian  86,  is  avail¬ 
able  from  Carnegie-Me 1  Ion  University  (Reference  21). 

Gaussian  86  is  a  connected  system  of  programs  for  performing  ab  initio 
molecular  orbital  (MO)  calculat ions .  It  represents  further  development  of  the 
Gaussian  70,  Gaussian  76,  Gaussian  80  and  Gaussian  82  systems  already 
published.  The  contributors  to  the  current  version  include:  M.  J.  Frisch,  J. 
S.  Binkley,  H.  B.  Schlegel,  K.  Raghavachar i ,  C.  F.  Melius,  R.  L.  Martin,  J.  J. 
P.  Stewart,  F.  W.  Bobrowicz,  C.  M.  Rohlfing,  L.  R.  Kahn,  D.  J.  Defrees,  R. 
Seeger,  R.  A.  Whiteside,  D.  J.  Fox,  E.  M.  Fleuder,  and  J.  A.  Pople.  Gaussian 
86  was  originally  implemented  on  the  chemistry  department  DEC  VAX  11/780 
computer  at  Carnegie-Mellon  University.  Since  then  this  program  has  been 
installed  on  a  number  of  different  computers. 

Gaussian  86  was  designed  with  a  transparent  input  data  stream,  making 
this  program  very  user  friendly.  All  of  the  standard  input  is  free-format  and 
mnemonic.  Reasonable  defaults  for  input  data  have  been  provided,  and  the 
output  is  intended  to  be  self-explanatory.  Mechanisms  are  available  for  the 
sophisticated  user  to  override  defaults  or  interface  their  own  code  to  the 
Gaussian  system.  In  this  respect,  we  intend  to  utilize  Gaussian  86  as  a 
fundamental  framework  for  several  applications.  Options  are  being  incor¬ 
porated  into  this  code  to  provide  capabilities  beyond  Hartree-Fock  and  various 
perturbation  theory  options. 

The  capabilities  of  the  Gaussian  86  system  include: 

a)  Calculation  of  one-  and  two-electron  integrals  over  s,  p,  d,  and  f 
contracted  gaussian  functions.  The  basis  functions  can  either  be 
cartesian  gaussians  or  pure  angular  momentum  functions  and  a  variety 
of  basis  sets  are  stored  in  the  program  and  can  be  requested  by 
name . 

b)  Self-consistent  field  calculations  for  restricted  closed-shell  (RHF), 
unrestricted  open-shell  (UHF),  and  open-shell  restricted  (ROHF) 
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Hartree-Fock  wave  functions  as  well  as  those  types  of  multiconf igura- 
tional  wavefunct ions  that  fall  within  the  Generalized  Valence  Bond- 
Perfect  Pairing  (GVB-PP)  formalism. 

c)  Evaluation  of  various  one-electron  properties  of  the  Hartree-Fock 
wavefunct ion ,  including  Mulliken  population  analysis,  multipole 
moments,  and  electrostatic  fields. 

d)  Automated  geometry  optimization  to  either  minima  or  saddle  points, 
and  numerical  differentiation  to  produce  force  constants,  polar¬ 
izabilities,  and  dipole  derivatives.  This  feature  can  be  used  to 
develop  minimum  energy  reaction  paths  along  a  complicated  many  dimen¬ 
sional  potential  energy  surface. 

e)  Correlation  energy  calculations  using  M<|>1  ler-Plesset  perturbation 
theory  carried  to  second,  third,  or  fourth  order. 

f)  Correlation  energy  calculations  using  configuration  interaction  (Cl) 
using  either  all  double  excitations  (CID)  or  all  single  and  double 
excitations  (CISD). 

g)  Correlation  energy  calculations  using  coupled  cluster  theory  with 
double  substitutions  (CCD). 

b)  Analytic  computation  of  the  nuclear  coordinate  gradient  of  the  RHF, 
UHF ,  ROHF ,  GVB-PP,  MP2 ,  CID  and  RCISD  energies. 

i)  Computation  of  force  constants  (nuclear  coordinate  second  deriva¬ 
tives),  polarizabilities,  hyperpolarizabilities,  dipole  derivatives, 
and  polarizability  derivatives  analytically  for  RHF  and  UHF  enegies, 
and  numerically  for  RHF,  UHF,  MP2,  CID  and  UCISD  energies. 

j)  Harmonic  vibrational  analysis. 

k)  Determination  of  intensities  for  vibrational  transitions  at  the  HF, 
MP2,  and  Cl  levels. 

l)  Testing  the  SCF  wave  functions  for  stability  under  release  of 
constraints . 

To  widen  the  appl icab il ity  of  Gaussian  86  to  open-shell  molecular 
species,  an  option  which  incorporates  a  spin-projected  variational  Cl  expan¬ 
sion  is  currently  under  development.  The  programs  for  constructing  the 
elements  of  the  hamiltonian  matrix  can  be  taken  from  our  standard  DIATOM  code. 
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Alternatively,  Cl  wave func t ions  can  be  constructed  using  the  unitary  group 
approach  presently  implemented  in  the  GAMESS  code. 

2.2  TRANSITION  PROBABILITIES. 

The  electronic  and  vibrational-rotational  wave func t ions  of  a  pair  of 

states  can  be  used  to  calculate  transition  probabilities.  If  two  molecular 

states  are  separated  in  energy  by  an  amount  Ae  =  hc^  (h  =  Planck's  constant, 

nm 

c  =  velocity  of  light,  v  =  frequency  in  wave  numbers),  the  semi-classical 
theory  of  radiation  (References  37  and  38)  yields  for  the  probability  of  a 
spontaneous  transition  from  an  upper  state  n  to  a  lower  state  m 


AE  n  m 


Here  A  is  the  Einstein  coefficient  for  spontaneous  transition  from  level  n 
nm 

m,  gn  is  the  total  degeneracy  factor  for  the  upper  state 


9n  =  (2-  8<r,A')(2S,+  |)(2J'  +  |) 


and  S  is  the  total  strength  of  a  component  line  in  a  specific  state  of 
nm 

polarization  and  propagated  in  a  fixed  direction.  A  related  quantity  is  the 
mean  radiative  lifetime  of  state  n  defined  by 
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the  summation  being  over  all  lower  levels  which  offer  allowed  connections. 
The  intensity  of  the  emitted  radiation  is 

1  nm  =  AEnm  Nn  Anm  (39' 


where  N  is  the  number  density  in  the  upper  state  n.  This  analysis  assumes 
n 

that-  all  degenerate  states  at  the  same  level  n  are  equally  populated,  which 

will  be  true  for  isotropic  excitation.  The  total  line  strength  S  can  be 

nm 


wntten  as  the  square  of  the  transition  moment  summed  over  all  degenerate 
components  of  the  molecular  states  n  and  m: 


=  I  M; 


where  j  and  i  refer  to  all  quantum  numbers  associated  collectively  with 
upper  and  lower  electronic  states,  respectively. 

In  the  Born-Oppenhe imer  approximation,  assuming  the  separability  of  elec 
tronic  and  nuclear  motion,  the  wavefunction  for  a  diatomic  molecule  can  be 
written  as 


^ vjmA  ^el  (R)v|/jm^(6,  x.  ^>) 


where  i}i  ,  (r,  R)  is  an  electronic  wavefunction  for  state  i  at  fixed  inter- 
nuclear  separation  R,  ^V(R)  is  a  vibrational  wavefunction  for  level  v  and 
0  (0,  x,  <}> )  refers  to  the  rotational  state  specified  by  electronic  angular 

JnA 

momentum  A,  total  angular  momentum  J  and  magnetic  quantum  number  M.  The 
representation  is  in  a  coordinate  system  related  to  a  space-fixed  system  by 
the  Eulerian  angles  (9,x,<{)).  The  transition  moment  M„  can  be  written,  using 
the  wavefunction  given  by  Equation  (41),  as 


M|i=/+ij'A'M' {m'  + «"}*(-/■  A“»"  dTedl»dTr 


The  subscripts  e,  v  and  r  refer  to  the  electronic,  vibrational  and  rota- 
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tional  wave funct ions ,  and  M  and  M  are  the  electronic  and  nuclear  electric 
dipole  moments,  respectively.  Integration  over  the  electronic  wavefunction, 
in  the  Born-Oppenhe imer  approximation,  causes  the  contribution  of  the  nuclear 


m 


moment  Mn  to  vanish  for  i  #  j.  The  electronic  dipole  moment  can  be  written 
(References  38  and  39)  in  the  form 


Me  =  -£er '  =- {  lerj.^s  (0lX.<*>)  (43) 

k  k 

where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordin¬ 
ates  r^  refer  to  a  molecule- fixed  system  and  (9,  X,  4>)  is  a  group  rotation 
tensor  whose  elements  are  the  direction  cosines  related  to  the  Eulerian  rota¬ 
tion  angles  (d.X.'M.  Using  bracket  notation,  Equations  (42)  and  (43)  can  be 
combined  to  yield  for  the  transition  moment 
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(44) 


The  matrix  elements  <J’  A’  M'j  ( 0 ,  x ,  i> )  |  J"  A"  M">  determine  the  group  se¬ 
lection  rules  for  an  allowed  transition  and  have  been  evaluated  for  many  types 
of  transitions  (References  40-42).  Summing  Equation  (44)  over  the  degenerate 
magnetic  quantum  numbers  M'  and  M" ,  we  have  from  Equation  (40) 
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where  j"/\"  is  the  Honl-London  factor  (References  43  and  44)  and 
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is  the  band  strength  for  the  transition.  Combining  Equations  (37),  (39)  and 
(45),  we  have  for  the  intensity  oi  a  single  emitting  line  from  upper  level  n: 
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where  Nj"  is  the  number  density  in  the  upper  rotational  state  J'  and 
“n  ■  (2-6o,A)  (2S '  +  1)  is  the  electronic  degeneracy.  Taking  an 
average  value  of  yMj,,  for  the  whole  band,  Equation  (48)  can  be  summed  to 
yield  the  total  intensity  in  the  (v1,  v")  band: 
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where  N^,  =  ^ ,  N^,  is  the  total  number  density  in  the  upper  vibrational  level 
v1  and  where  we  make  use  of  the  group  summation  property 
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Comparing  Equations  (39)  and  (48),  we  have  for  the  Einstein  spontaneous 
transition  coefficient  of  the  band  (v‘,  v") 
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Similarly,  the  lifetime  of  an  upper  vibrational  level  v1  of  state  n  can  be 


where  the  summation  runs  over  all  v"  for  each  lower  state  m.  Equation  (50) 
can  be  cast  in  the  computational  form 
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where  AE°  V..  and  pn  V,,  are  in  atomic  units.  It  is  also  often  convenient  to 
m  v  m  v 

relate  the  transition  probability  to  the  number  of  dispersion  electrons  needed 
to  explain  the  emission  strength  classically.  This  number,  the  f-number  or 
oscillator  strength  for  emission,  is  given  by 
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Tlie  inverse  process  of  absorption  is  related  to  the  above  development 
through  the  Einstein  B  coefficient.  Corresponding  to  Equation  (39),  we  have 
for  a  single  line  in  absorption 
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where  K(v)  is  the  absorption  coefficient  of  a  beam  of  photons  of  frequency 


v  and 
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is  the  Einstein  absorption  coefficient  for  a  s ingle  line.  Summing  over  al 
1  ines  in  the  band  (v",  v'),  assuming,  an  average  band  frequency,  we  obtain 
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where  N  „  =  N ...  is  Che  total  number  density  in  the  lower  vibrational  state 
v  J  J 

v" .  Corresponding  to  Equations  (52)  and  (53)  we  can  define  an  f-number  or 
oscillator  strength  for  absorption  as 


f 


mn,v'V 


—  nv' 
2^  AEmv" 

3Ti2e2ujm 


nv 

Pmv" 


(57) 


In  computational  form,  Equation  (57)  becomes 
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where  &E  ,,  and  p  ,,  are  in  atomic  units.  Combining  Equations  (50)  and 
m  v  rm  v 

(53)  and  comparing  with  Equation  (58),  we  see  that  the  absorption  and  emission 
f-numbers  are  related  by 
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Some  caution  must  be  observed  in  the  use  of  f-numbers  given  either  by  Equa¬ 
tion  (53)  or  (57)  since  both  band  f-numbers  and  system  f-numbers  are  defined 
in  the  literature.  The  confusion  arises  from  the  several  possible  band  aver¬ 
aging  schemes  that  cap.  be  identified. 

An  integrated  abs  rption  coefficient  (density  corrected)  can  be  defined 
from  Equation  (56)  as 
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where  the  exponential  factor  corrects  for  stimulated  emission.  Equation  60 
can  be  written  in  terms  of  the  absorpt ion  f-number  as 
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Using  hc/k  -  1.43880  cm-K° ,  we  obtain  a  computational  formula  for  the 
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integrated  absorption  coefficient  as  Svm  vi  (cm  atm  )  = 
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The  total  integrated  absorption  is  found  from 
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where,  under  normal  temperature  conditions,  only  the  first  few  fundamental s 
and  overtones  contribute  to  the  summations. 

The  developments  given  above  are  rigorous  for  band  systems  where  an  aver 
age  band  frequency  can  be  meaningfully  defined.  Further  approximations,  how¬ 
ever,  are  often  made.  For  example,  the  electronic  component  of  the  dipole 


transition  moment  can  be  defined  as 


This  quantity  is  often  a  slowly  varying  function  of  R  and  an  average  value 
can  sometimes  be  chosen.  Equation  (46)  can  then  be  written  approximately  in 
factored  form  as 
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where  q  .  .. 
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the  square  of  the  vibrational  overlap  integral,  is  called  the 


Franck-Condon  factor.  D..  is  evaluated  at  some  mean  value  of  the  internuclear 
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separation  R.  In  addition,  it  is  sometimes  possible  to  account  for  a  weak  R- 
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dependence  in  M  by  a  Taylor  series  expansion  of  this  quantity  about  some 
reference  value,  usually  referred  to  the  (0,  0)  band.  We  have 
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Substituting  into  Equation  (66)  and  integrating  yields 
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where 


( R  v '  v "  -  Ra/3)  = 


<V|lR-Ra;S)|v'> 

<V'|  V"> 


(68) 


is  the  R-centroid  for  the  transition  and 


is  the  R^-centroid.  Note  that  this  last  term  differs  (to  second  order)  from 
the  square  of  the  R-centroid.  An  alternate  procedure  can  be  developed  by 
evaluating  Equation  (64)  at  each  R-centroid,  Rvtvn.  Then 
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Equation  (70)  assumes  that  the  vibrational  wavefunction  product  0  41 

v 

behaves  like  a  delta  function  upon  integration, 


W ^v,=  S(R~  Rv'v")  <Y'|v'^> 


(71) 


The  range  of  validity  of  Equation  (70)  is  therefore  questionable,  particu¬ 
larly  for  band  systems  with  bad  overlap  conditions  such  as  oxygen  Schumann- 
Runge.  The  range  of  validity  of  the  R-centroid  approximation  has  been  exam¬ 
ined  by  Fraser  (Reference  45). 

The  final  step  in  calculating  transition  probabilities  is  the  determina¬ 
tion  of  $j^(R),  the  electronic  dipole  transition  moment,  for  the  entire  range 
of  internuclear  separations,  R,  reached  in  the  vibrational  levels  to  be  con¬ 
sidered.  This  can  be  expressed  in  terms  of  the  expansion  of  Equation  (14) 
as 


i  <R)  =Z  cj»cj,  <*m(R)|  MeL„(R)> 
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where  cuJ  and  cj  are  coefficients  for  and  ^  ,  respectively. 

An  analysis  similar  to  that  yielding  Equation  (10)  and  (11)  gives 
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The  spatial  integral  in  Equation  (73)  reduces  to  one-electron  integrals 

equivalent  to  overlap  integrals,  and  the  evaluation  of  Equation  (73)  can  be 

carried  out  by  the  same  computer  programs  used  for  Equation  (11).  Programs 

for  evaluating  $..( R)  in  Equation  (72)  have  been  developed  at  UTRC  and  exam- 
Ji 

pies  of  their  application  have  appeared  in  the  literature  (Reference  8). 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have 
a  strong  R-dependence  and  R-centroid  or  other  approximations  will  be  invalid. 
A  direct  evaluation  of  Equation  (46)  would  therefore  be  required  using  the 
fully-coupled  system  of  electronic  and  vibrational  wavefunct ions  to  properly 
account  for  the  source  of  the  band  perturbations. 
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SECTION  3 


DISCUSSION  OF  RELATIVISTIC  METHODS 


For  heavy  atoms  (Z  >  30),  and  molecular  systems  built  from  heavy  atoms, 

relativistic  effects  become  increasingly  important  and  should  be  taken  into 

account  in  the  calculation  of  the  radial  wavefunct ions .  The  implementation  of 

relativistic  effects  into  atomic  and  molecular  computer  codes  is  only  fairly 

recent  owing  to  the  increased  complexities  introduced  in  the  self-consistent 

field  (SCF)  procedure  and  the  greatly  increased  computer  time  required  for 

such  calculations.  Compared  with  the  non-relat iv ist ic  case,  the  Dirac- 

Hartree-Fock  (DHF)  method  requires  that  two  radial  functions,  G  , corre- 

nl  j 

sponding  to  the  large  component  and  F  .  corresponding  to  the  small  component 

n  1  j  , 

must  be  calculated  for  each  of  the  two  possible  j  values.  Thus,  the  numerical 
work  of  a  DHF  relativistic  treatment  is  increased  by  nearly  a  factor  of  four 
over  the  nonrelat ivist ic  case,  exclusive  of  increased  complexities  in  evalua¬ 
tion  of  the  terms  of  the  Hamiltonian.  In  view  of  this,  methods  that  have  been 
developed  to  date  for  molecular  systems  have  involved  the  use  of  model  poten¬ 
tials  to  represent  relativistic  effects. 

In  the  calculation  of  the  internal  energy  of  a  molecular  system  comprised 
of  n  electrons  and  N  nuclei,  and  considering  only  electrostatic  interactions 
between  the  particles,  we  have  for  the  total  Hamiltonian 
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where 


where  m  ,  m  ,  M  ,  are  the  masses  of  the  electron,  atom  a  and  combined  system 
e  a  T 

mass,  respectively.  Now  since  the  ratios  m  /m  and  m  /M„,  are  both  small, 
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(2xl0-c>  -  5xl0-1+)  we  can  effect  a  separation  of  the  electronic  and  nuclear 

coordinates  treating  the  total  wavefunction  as  a  product  of  a  nuclear  and  an 

electronic  part.  We  have 
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where  if  (r  ,  R  )  is  an  electronic  wavefunction  parametric  in  the  nuclear 
k  ~n  N 

coordinates  as  given  in  Equation  (76)  and  are  nuclear  motion  wavefunc- 

t  ions  which  satisfy  (neglecting  terms  of  the  order  of  mg/rt^) 
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The  cross  terra  in  7^  *  can  be  eliminated  by  a  proper  change  of  varia¬ 
bles  and  Equation  (77)  then  reduces  to  a  3N-3  dimensional  Schrodinger 
equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow 

relative  to  the  electron  velocity,  this  decoupling  of  electronic  and  nuclear 

motion  is  valid  and  is  referred  to  as  the  adiabatic  approximation.  Equation 

(76)  thus  defines  an  electronic  eigenstate  ii>,  (r  ,  R  ) ,  parametric  in  the 

k  ~n  ~N 

nuclear  coordinates,  and  a  corresponding  eigenvalue  E^(R^)  which  is  taken  to 
represent  the  potential  energy  curve  or  surface  corresponding  to  state  k. 

In  the  usual  ab  initio  method  for  calculating  the  electronic  properties 
of  a  molecular  system,  one  starts  from  a  zero-order  Hamiltonian  that  is  exact 
except  for  relativistic  and  magnetic  effects,  and  which  involves  the  evalua¬ 
tion  of  electronic  energies  and  other  relevant  quantities  for  wave func t ions 
that  are  properly  antisymmetrized  in  the  coordinates  of  all  the  electrons. 
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For  a  system  containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian 
depends  paramet r ical ly  on  the  nuclear  positions  and  is  of  the  form 
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where  z.  and  R.  are  the  charge  and  position  of  nucleus  i,  r.  is  the  position 
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of  electron  j,  and  7^  is  the  Laplacian  operator  for  electron  j.  All  quanti¬ 
ties  are  in  atomic  units,  i.e.  lengths  in  bohrs,  energies  in  hartrees 
(1  hartree  =  2  Rydbergs). 

In  addition  to  the  electrostatic  contribut  ion,  tTf-  ,  the  complete  Hamilton¬ 
ian  should  contain  additional  terms  which  correct  for  magnetic  interactions 
and  relativisitic  effects.  These  correction  terms  may  be  of  importance  in 
several  applications.  These  include: 

(1)  calculation  of  the  probability  of  making  a  transition  from  one  quan¬ 
tum  state  to  another  in  high-momentum  collisions  such  as  those  that 
can  occur  in  hot  atom  or  heavy  atom  chemical  dynamics  experi¬ 
ments  ; 

(2)  determination  of  the  interaction  energy  in  heavy  nuclei  systems  such 
as  Cs2  and  LiO+,  which  exhibit  open-shell  structure  on  both  nuclei  at 
infinite  internuclear  separations; 

(3)  calculation  of  the  intermolecular  forces  between  free  radicals, 
electronically  excited  states  of  molecules  with  open-shell  struc¬ 
ture,  and  long  molecular  conformations  of  possible  biological 
interest . 

3.1  BRE IT-PAULI  HAMILTONIAN. 

The  relativistic  correction  terms  to  the  usual  electrostatic  Hamiltonian 
have  been  derived  through  order  ,  where  a  is  the  fine  structure  constant, 
and  are  often  referred  to  as  the  Breit-Pauli  (Reference  46)  Hamiltonian  terms. 
This  Hamiltonian  has  been  derived  by  Bethe  and  Salpeter  (Reference  47)  for  a 
two-electron  system  and  has  been  generalized  to  the  many-e lec t ron  system  by 
Hirschfelder ,  et  al  (Reference  48)  and  Itoh  (Reference  49).  In  the  absence  of 
external  electric  or  magnetic  fields  we  can  represent  these  correction  terms 
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as  follows .  Let  s . 

J 

moment  of  electron 
tonian,  correct  to 


and  p.  =  —  A  .  denote  the  operators  for 
J  i  J 

j,  respectively.  Then  the  generalized 
terms  of  Q(<*^ /M)  ,  can  be  written  as: 


the  spin  and  linear 
Breit-Pauli  Hamil- 


-*BP=  A+Al  +  As  *  As  +  A  +  A 


(79) 


where^fg  is  given  by  Equation  (78)  and  the  correction  terms  can  be 
expressed  as  follows: 


Jl  = — 4"  Y  1  t  fr2  "p.  .  p  +  r  •  (r  •  p)  p  1 
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The  tirst  correction  term, ^  ,  represents  the  magnetic  orbit-orbit  coupling 

Li  Li 

terms  of  the  electrons  arising  from  the  interaction  of  the  magnetic  fields 
created  by  their  motion.  The  second  term,,^^,  gives  the  spin-spin  magnetic 
coupling  terras  which  are  often  quite  appreciable.  For  r  =0,  only  the 
delta-function  contribution  survives  which  represents  the  Fermi-contact  spin 
interaction.  The  third  term,o^  ,  is  usually  the  largest  in  magnitude  and 

Ld 

represents  the  spin-orbit  interaction  between  the  spin  and  magnetic  moment  of 

each  electron  and  the  spin-other  orbit  interaction,  which  represents  the 

coupling  of  the  spin  of  one  electron  with  the  magnetic  moment  of  a  different 

electron.  The  term, corrects  for  variation  of  the  electron  mass  with 

velocity  and  the  term)C>^Q,  represents  electron  spin  terms  identified  by 

Dirac  which  appear  to  have  no  classical  analogue. 

Asido  from  the  spin-orbit  term,  ^  ,  usually  only  the  last  term,  , 

Li  S  0 

(often  called  the  Darwin  correction  term)  and^p,  the  mass-velocity  term, 
are  retained  in  the  Hamiltonian,  yielding  the  so-called  Pauli  approximation 
(Reference  47) . 

The  eigenfunctions  of  the  Hamiltonian  represented  by  Equation  (79)  are 
four-component  Dirac  spinors  which  may  be  expressed  as: 


Pnk^ 

'<Wr! 


where  X,  (9>  4> )  are  products  of  spherical  harmonics  and  Pauli  spinors  and 
(on 

P  ^(r) ,  Qnk(r)  represent,  respectively,  the  large  and  small  components  of  the 
radial  wave  equation.  The  exact  solution  of  the  Breit-Pauli  Hamiltonian  has 
only  been  given  for  one-  and  two-electron  atomic  systems  (Reference  47)  owing 
to  the  complexity  of  the  operators  for  the  general  n-electron  case.  For  a 
molecular  system,  Kolos  and  Wolniewicz  (Reference  50)  have  calculated  the 
relativistic  corrections  to  H^  using  Equation  (79)  to  (P(a^).  No  heavier  mole 
cular  systems  have  been  treated  using  the  full  Breit-Pauli  Hamiltonian. 

3.2  APPROXIMATE  TREATMENTS. 

Although  the  Breit-Pauli  Hamiltonian  given  in  Equation  (79)  can  formally 
be  employed  in  a  molecular  system,  both  the  multiplicity  of  terms  and  the 


difficulty  of  evaluation  of  the  resultant  molecular  integrals  has  precluded 
its  general  use  to  date.  For  atomic  systems,  various  approximate  methods  of 
solution,  within  a  Hartree-Fock  or  multiconfiguration  Hartree-Fock  framework, 
have  been  proposed  for  atoms  (References  51-55).  In  most  of  these  methods,  a 
restricted  Hamiltonian  which  includes  only  the  one-electron  Dirac  terms  is 
usually  employed.  The  contributions  of  the  Breit  operators  for  spin-magnetic 
interactions  and  velocity  retardation  are  then  calculated  as  first-order  per¬ 
turbations  using  the  zeroth-order  Dirac  relativistic  wavefunctions . 

An  even  more  approximate  method  for  incorporating  the  major  relativistic 
effects  has  been  proposed  by  Cowan  and  Griffin  (Reference  56).  In  this 
method,  the  mass-veloc ity ( cMp)  and  Darwin (^q)  terms,  written  in  terms  of 
the  Pauli  equation  for  one-electron  atoms,  are  simply  added  to  the  usual  non- 
relativistic  Hamiltonian  operator.  In  addition,  the  spin-orbit  terms,  , 
are  omitted,  thereby  reducing  the  system  of  equations  to  a  single  form  repre¬ 
senting  the  description  of  the  major  component  wavefunction,  P  (r) ,  evaluated 
at  the  center-of-gravity  of  the  spin-orbit  states.  The  rationale  for  this 
approximation  lies  in  the  observation  that  detailed  atomic  calculations  using 
the  complete  DHF  method  have  indicated  that,  even  for  an  atom  as  heavy  as 
uranium,  less  than  1  percent  of  the  total  charge  is  described  by  the  small 
component  radial  wavefunctions. 

The  resulting  equations  have  the  form: 


l  i  (*i+l) 


+  V1  (r)  +  (r)  +  hd'  (r)  G'nl  (r)  =  «'ni  G*n£(r) 


(86) 


where  the  mass-velocity  and  Darwin  terms  take  the  form 


H'm  (r) 


(87) 


(88) 
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and  a  ~  1/137.036  is  the  fine  structure  constant.  The  spin-orbit  term  is 
omitted  in  Equation  (86)  aiui  l'iius  tiie^n  equation.-,  represent  center  of  gravily 
radial  functions  averaged  over  the  two  possible  total  angular  momentum  quantum 
numbers.  Equation  (86)  represents  (apart  from  the  neglect  of  spin-orbit 
effects)  the  relativistic  corrections  to  first  order  in  «  .  A  more  accurate 
analysis  of  heavy  atom  energy  levels  and  spectra  is  available  through  the  use 
of  the  radial  functions,  G1  ^(r),  found  from  Equation  (86),  and  a  first-order 
perturbation  calculation.  Cowan  and  Griffin  (Reference  56)  have  illustrated 
the  utility  and  accuracy  of  such  an  approach. 

Recently,  Wood  and  Boring  (Reference  57)  have  adapted  this  approximate 
relativistic  method  to  the  local  exchange  problem  and  have  implemented  the 
solution  of  Equation  (86)  within  the  context  of  the  multiple  scattering  X« 
method  (Reference  26).  The  central  field  Hamiltonian  is  modified  to  include 
mass-velocity  and  Darwin  terms,  given  by  Equations  (87)  and  (88),  in  the 
sphere  surrounding  each  atomic  center.  The  intersphere  region  in  the  multiple 
scattering  approach  (constant  potential  region)  is  treated  nonrelat ivist  ical ly 
since  charge  in  this  region  is  far  from  a  nucleus  and  is  screened  by  the 
charge  concentrated  around  the  atomic  centers.  The  matching  conditions  for 
continuity  of  the  wavefunction  at  the  sphere  boundaries  permits  any  necessary 
charge  transfer  between  the  relativistic  intra-atomic  regions  and  the  nonrela- 
tivistic  interatomic  constant  potential  regions.  For  an  atom,  the  Wood-Boring 
treatment  reduces  to  the  Dirac-Slater  local  exchange  method,  but  with  the 
neglect  of  spin-orbit  terms. 

The  implementation  of  Equation  (86)  into  existing  nonrelat  ivist  ic  multi¬ 
ple  scattering  molecular  codes  is  facilitated  by  a  change  in  the  dependent 
variable,  Gl^(r),  to  eliminate  the  first  derivative  of  the  wavefunction, 
illustrated  in  Equation  (88).  Hie  usual  Numerov  method  of  solution  can  then 
be  applied  to  the  central  field  problem;  the  only  new  requirement  being  the 
numerical  tabulation  of  the  first  and  second  derivatives  of  the  potential  at 
each  grid  point  in  the  integrations.  These  derivatives  are  computed  only  once 
for  each  complete  SCK  cycle  and  thus  the  total  required  computer  time  for  a 
typical  problem  is  not  significantly  increased  as  compared  with  a  nonrelat iv- 
istic  calculation.  A  complete  self-consistent  program  incorporating  this 
method  has  been  developed  at  UTRC .  Our  code  has  been  tested  by  repeating 
calculations  for  the  U  and  I’u  atoms  (Reference  56),  where  we  find  excell.:. it 
agreement  with  the  more  exact,  but  cumbersome,  Dirac-Slater  calculations. 


Results  for  molecular  calculations  have  recently  been  reported  by  Boring  and 
Wood  for  UF k  and  (References  58  and  59).  These  calculations  were  car¬ 

ried  out  to  illustrate  the  shifts  in  the  valence  levels  for  such  systems  re¬ 
sulting  from  relativistic  effects.  The  total  energy  was  not  of  principal 
concern . 

We  have  recently  reported  (Reference  60)  the  first  all-electron  calcula¬ 
tion  of  the  potential  energy  curves  for  a  molecule  (Hg^+)  built  from  atoms 
which  exhibit  significant  relativistic  effects.  This  study  illustrated  that 
reliable  total  energies  are  obtainable  through  a  relativistic  multiple  scat¬ 
tering  density  functional  treatment,  provided  care  is  taken  to  optimize  poten¬ 
tial  match  and  overlap  criteria  for  such  systems.  This  study  formed  the  basis 
of  the  computational  scheme  that  we  have  employed  here  for  the  uranium/oxygen 
system. 

3.3  EFFECTIVE  CORE  MODELS. 

It  is  well  known  from  chemical  experience  that  the  outermost  valence 
electrons  contribute  most  to  determining  the  chemical  properties,  especially 
spectroscopic  properties,  of  molecules.  The  core  electrons  remain  essentially 
unchanged  from  their  atomic  form  except  for  internuclear  separations  of  the 
order  of  the  charge  radii  of  the  outer  core  region  or  less,  wherein  core 
polarization  effects  may  become  important.  Since  the  computational  time 
required  for  ab^  initio  calculations  of  electronic  structure  goes  up  at  least 
quadratically  with  the  number  of  electrons  in  the  system,  there  have  been  many 
attempts  to  replace  the  more  tightly  bound  core  electrons  with  simple  one- 
electron  effective  potentials  (References  61-69).  Concurrent  with  elimination 
of  an  explicit  treatment  of  the  core  electrons,  a  transformation  of  the 
valence  orbital  basis  is  required  to  insure  that  the  lowest  valence  orbital  of 
each  symmetry  has  a  nodeless  radial  form,  since  it  is  well  known  that  the 
lowest  energy  eigenfunction  for  a  local  potential  must  be  nodeless  (Reference 
70). 

Typical  of  the  several  effective  core  models  that  have  been  reported  is 
that  due  to  Kahn,  et  al  (Reference  68)  whereby  an  effective  core  potential  is 
described  in  terms  of  angularly  dependent  projection  operators  as 

Uco,e  =  UL“"(r|  +  I  I  |im>[u,cor'(r)-uLcor'(r)]<|m|  (89) 


where  L  is  taken  at  least  as  large  as  the  highest  angular  momentum  orbital 
occupied  in  core.  The  term  U  C°re(r)  represents  the  effective  Coulomb  and 

lj 

exchange  potential  felt  by  the  valence  electrons.  The  second  term  essentially 
accounts  for  the  repulsive  potential  between  valence  and  core  electrons  for 
each  symmetry  £.  The  only  non-local  character  exhibited  by  a  potential  of  the 
form  of  Equation  (89)  arises  from  the  H -dependence  which  '■an  be  cast  in  terms 
of  one-electron  integrals  between  the  core  and  valence  orbitals.  Explicit 
two-electron  terms  connecting  core  and  valence  orbitals  are  thus  avoided  which 
greatly  simplifies  the  calculation  of  matrix  elements  of  the  effective  Hamil¬ 
tonian.  The  potential  given  by  Equation  (89)  can  be  compared  with  the  gener¬ 
alized  Ph i l 1 ips-Kl e inman  pseudo-potential  (Reference  63). 

UCore  =  I  (2Jr  -  Kc)  +  VC0  (90) 

c  c  c 

where  Jc  and  Kc  represent  the  core  orbital  Coulomb  and  exchange  operators 
CO 

and  V“  is  a  complicated  non-local  operator  which  guarantees  core-valence 

orthogonality.  Since  the  P-K  core  orbitals  must  simultaneously  be  eigenfunc- 

CO 

tions  of  both  the  core  and  valence  Hartree-Eock  Hamiltonians,  V  ,  in  general, 
contains  complicated  two-electron  terms  and  limits  the  usefulness  of  Equation 
(90)  over  a  full  ab  initio  treatment. 

The  prescription  of  Kahn  can  be  implemented  by  analytically  fitting  a 
node  less  pseudo-orb i ta  1 ,  x  to  a  linear  combination  of  numerical  or  analytic 
Hartree-Fock  orbitals  determined  from  a  full  self-consistent  treatment  of  the 
core  electrons,  using,  for  example,  tne  multiconfiguration  Hartree-Fock  code 
of  Froese-Fischer  (Reference  71).  The  components,  UfC°ri  (r),  of  Equation 
(89)  are  then  defined  implicitly  from  the  Schrodinger  equation 


[-  y  --f  +  U,COre(r)  +  2JvQ1-Kvq1]  Xnf  =tni  Xni  (<M  ) 
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Equation  (91)  can  be  extended  to  relativistic  systems  in  several  ways. 
Kahn,  et  al  (Reference  72)  suggest  an  approximate  treatment  of  adding  only 
the  mass-velocity  and  Darwin  terms  to  the  usual  electrostatic  Hamiltonian  and 
to  determine  approximate  HF  orbitals  in  the  manner  prescribed  by  Cowan  and 
Griffin  (Reference  56).  Equation  (91)  is  then  used  to  determine  an  effective 
U««  (r)  such  that  e  ^  are  the  eigenvalues  of  the  CG  approximate  relativis¬ 
tic  solution  and  X  „  are  curve-fitted  to  the  CG-HF  orbitals.  In  this  treat- 

nx 

ment ,  the  Xn^  represent  approximate  solutions  to  the  major  component  wave  func¬ 
tion,  Pn£ ,  determined  at  the  center-of-gravity  of  the  spin-orbit  states. 

Lee,  et  al  (Reference  69)  adopt  a  somewhat  more  complicated  treatment  in 
which  the  spin-orbit  operator  is  added  to  the  usual  electrostatic  Hamiltonian, 
in  addition  to  the  mass-velocity  and  Darwin  terms  retained  in  the  Cowan- 
Griffin  treatment.  The  large  component  eigenfunctions  of  a  full  Dirac- 
Hartree-Fock  treatment  of  the  atom,  as  given,  for  example,  by  Desclaux  (Refer¬ 
ence  73)  are  then  curve-fitted  in  a  manner  similar  to  the  Kahn  treatment  but 

include  the  additional  index  for  the  particular  spin-orbit  state,  X  ...  Use 

nx  j 

of  these  eigenfunctions  in  a  molecular  system  fits  more  naturally  into  a  (J-J) 
coupling  scheme  whereas  the  X  £  determined  using  the  Kahn  method  are  more 
easily  represented  using  A-s  coupling. 

Although  these  effective  core  models  can  often  accurately  describe  an 
atomic  eigenvalue  sequence,  including  even  high-lying  electronically  excited 
states  (Reference  74),  there  are  inherent  difficulties  in  their  application  to 
molecular  environments,  where  the  maximum  angular  momentum  component  of  the 
valence  shell  orbitals  may  often  exceed  the  highest  l -value  component  retained 
in  Equation  (92).  This  is  particularly  true  for  valence  orbitals  which  exhi¬ 
bit  strong  changes  from  atomic  form  through  hybridization  with  higher  angular 
momentum  orbitals  or  through  the  addition  of  more  compact  polarization  terms. 
In  either  case,  since  the  relativistic  terms  are  now  all  buried  in  a  fixed 
rather  than  a  dynamic  relativistic  operator,  only  static  core  effects  are 


imposed  in  determining  the  shape  of  the  valence  molecular  orbitals.  Relativ¬ 
istic  effects  between  valence  electrons  and  shielding  effects  of  the  core  by 
the  valence  electrons  are  therefore  neglected  in  these  effective  core  treat¬ 
ments.  In  addition,  the  models  obviously  break  down  completely  when  the 
nuclei  are  brought  together  to  dimensions  such  that  core  overlap  and  polariza¬ 
tion  effects  become  significant.  Unfortunately,  calculations  to  date  seem  to 
indicate  that  such  effects  begin  to  set  in  for  internuclear  separations  of  the 
order  of  equilibrium  bond  lengths. 


states  will  be  optically  connected  and  many  pre-dissociation  paths  of  the 
type: 

U0'<+  +  hv  -►  IK)’*"  +  U  +  o  (94) 

*  +  . 

are  possible.  Here  UO  i'  a  v  ibrat Lonal 1 y  excited  low-mu  1 1  lplet  state  of  UO 
and  UO  is  a  dissociating  state.  The  predicted  optical  absorption  should  be 
strong  since  the  transfer  is  from  highly  ionic  states  to  neutral  valence 
states  of  UO.  Since  UO  has  a  large  dissociation  energy  (7.87  eV),  both  one 
photon  and  two  photon  solar  excitation  processes  are  possible. 

A  vibrational  analysis  of  the  ft  =  5  ground  state  of  UO  was  carried  out 
using  a  Hulbert-Hirshfelder  (Reference  77)  fit  to  our  calculated  potential 
curves.  This  fit  yields  an  equilibrium  internuclear  distance  of  1.89  A  and  a 
fundamental  vibrational  constant  (ug)  of  859  cm-1.  The  spin-orbit  interac¬ 
tion  was  calculated  using  U+2  atomic  splitting  parameters.  No  explicit  two- 
center  effects  are  included.  Our  calculated  spectroscopic  data  are  compared 
in  Table  1  with  the  theoretical  work  of  Krauss  and  Stevens  (Reference  78),  the 
recent  rotationally  resolved  experimental  studies  of  Heaven  and  Nicolai 
(Reference  79),  the  recent  rotational  analysis  reported  by  Kaledin,  et  al 
(Reference  80),  and  estimates  based  on  experimental  data  for  similar  systems. 
The  agreement  is  well  within  the  uncertainty  of  the  calculations  or  experi¬ 
mental  estimates.  The  theoretical  studies  predict  a  5I5  or  5I4  ground  state 
for  UO  whereas  the  experimental  data  of  both  Heaven  and  Nicolai  and  Kaledin, 
et  al  suggest  a  5I4  ground  state.  Recent  ligand  field  model  calculations  by 
Dulick  (Reference  81)  predict  a  ground  state  for  UO,  in  substantial 

disagreement  with  previous  theory  and  experiments.  The  crystal  field  model 
totally  neglects  overlap  and  back  charge  transfer  effects  in  the  U+20-2  model, 
both  of  which  have  been  found  to  be  of  importance  in  more  exact  studies 
carried  out  within  a  full  molecular  framework.  The  character  of  all  of  the 
low-lying  multiplets  of  UO  is  similar  however,  and  this  apparent  discrepency 
should  not  affect  our  conclusions  about  optical  or  LWIR  absorptions. 

An  analysis  of  the  LWIR  emission  from  UO  was  carried  out  based  on  a 
ground  ^1  electronic  state.  These  calculations  should  also  be  representative 
of  the  LWIR  emission  from  other  low-lying  electronic  states  since  they  exhibit 
similar  ionicities.  Our  calculated  f-numbers,  including  fundamentals  and 
overtones,  were  reported  in  DNA-TR-82-1 59  (Reference  75).  These  data 
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(v1  >  v")  were  given  for  the  lowest  30  vibrational  levels.  For  UO,  our 

calculated  f-number  for  the  1-0  transition  is  4.86  x  10-3  at  A  =  12.04  u. 

4.2  U0+. 

Detailed  searches  of  several  symmetries  of  U0+  were  carried  out  to  deter¬ 
mine  the  ground  molecular  state  of  this  system.  Our  calculations  indicate 
that  the  lowest  symmetry  of  U0+  is  derived  from  the  (4,  S)  coupled  4I  state 
and  has  the  following  principal  molecular  orbital  occupancy: 

^1  [la2  2a2  3a2  4a2  lw^  2 ^  !  3u  l<j>  15  ;J  (95) 

t  I 

Again  we  found  that  the  J  3tt  1£  16  J  group  is  quartet  coupled  in  the  ground 
state  but  a  second  manifold  of  states  for  the  U0+  system,  which  exhibit  doub¬ 
let  coupling  of  these  electrons,  lies  about  2  eV  above  the  ground  state.  An 
apparent  gap  in  the  density  of  states  for  U0+  is  found  between  these  two 
groups.  Calculations  for  UO+  proved  to  be  much  more  complex  than  UO  owing  to 
the  presence  of  at  least  two  low-lying  dissociation  limits  of  U+  +0.  In  a 
molecular  framework,  the  ionic  U++0~  and  U+;50~2  structures  will  mix  in  all 
multiplicities  with  molecular  ion  states  arising  from  U+[HIJ  +  0 (_ 3 P ]  ,  U+[bLj  + 
0[JpJ,  and  +  0  [  3  P  J  . 

A  vibrational  analysis  of  the  ft  =  9/2  ground  state  of  U0+  was  carried  out 
using  a  Hulbert-Hirschfelder  (Reference  77)  fit  to  our  calculated  potential 
curves.  The  spin-orbit  splittings  were  derived  from  atomic  parameters  for  the 
U+;J  ion.  This  fit  yields  an  equilibrium  internuclear  distance  of  1.84  A  and  a 
fundamental  vibrational  constant  of  890  cm-1.  These  data  are  compared  in 
Table  1  with  other  calculated  estimates,  since  there  are  no  experimental  data 
available.  The  agreement  between  our  work  and  that  of  Krauss  and  Stevens 
(Reference  78)  is  less  satisfactory  than  in  the  case  of  UO  but  still  well 
within  the  uncertainty  of  the  several  calculations. 

A  perturbative  treatment  for  calculating  the  density  of  states  in  uranium 
molecules  is  available  through  the  use  of  ligand  field  theory.  Recent  studies 
by  Dulick  predict  that  the  ground  state  of  U0+  has  ft  =  3/2,  in  contrast  to  the 
effective  core  hamiltonian  calculations  of  Krauss  and  Stevens  (Reference  78) 
and  our  relativistic  density  functional  calculations,  both  of  which  predict 
that  the  ground  state  has  ft  =  9/2.  Again,  the  neglect  of  important  molecular 
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effects  in  the  ligand  field  model  apparently  results  in  a  bias  toward  lower  fj 
values.  Our  calculated  spin-orbit  splitting  for  the  5f  electrons  in  U+ 3  u+5 

is  given  in  Table  2  where  a  comparison  with  earlier  studies  by  Wood  and  Boring 
(Reference  57)  is  given.  The  overall  agreement  for  these  atomic  ions  is 
excellent.  We  have  carried  out  calculations  of  the  excitation  energies  of 
several  configurations  of  U+3.  These  data  clearly  indicate  that  the  4I 
component  of  5f3  is  low-lying  and  that  the  lowest  doublet  manifold  lies  ~  2.2 
eV  above  the  ground  state.  If  this  splitting  carries  over  to  the  U0+  ion 
without  much  change,  the  following  solar  pumped  process  is  possible: 

U0+[4l]  +  hv  *  UO+*[2H]  (96) 

e  +  U0+*  [2H  ]  -*•  U  [5L  ]  +  0 [3P  ]  (97 ) 

Dissociative-recombination  is  not  energetically  possible  from  the  ground  4I 
state  of  U0+  but  it  is  possible  from  the  excited  and  higher  states.  The 
optical  connections  between  the  manifolds  of  U0+  states  will  be  studied  in 
future  work. 

An  analysis  of  the  emission  characteristics  for  the  ground  state  of  U0+ 
indicates  an  oscillator  strength  for  emission  (f  1  Q)  of  5.17  x  10" 5  at  X  = 

11.3  y.  A  complete  analysis  of  our  calculated  LWIR  emission  for  UO+  was  given 
in  DNA  TR-82-159  (Reference  75).  Our  calculated  LWIR  emission  for  U0+  is 
typical  of  that  for  a  highly  ionic  metal  oxide.  We  predict  strong  emission 
from  the  fundamentals  of  U0+  in  the  wavelength  region  11  ~  14  y.  Since  this 
system  exhibits  weak  anharmonic ity ,  we  find  the  overtones  down  in  intensity  by 
several  orders  of  magnitude.  However,  the  first  excited  state  of  U0+  ( 4H) 
lies  at  ~  2800  cm-1  in  our  calculations  with  a  predicted  electronic  oscillator 
strength  of  ~  1  x  10~5  for  the  4I  -  4H  transition.  The  electronic  and  vibra¬ 
tional  manifolds  for  U0+  are  highly  overlapped  above  the  second  vibrational 


level  of  the  ground  5Ig/2  state.  Since  the  density  of  electronic  states  of 
U0+  is  very  large  above  ~  2.2  eV,  we  predict  that  strong  solar  pumping, 
followed  by  intense  radiation  in  the  region  0.6  <  X  <  11.3  y  should  occur  fcr 
this  system.  This  conclusion  is  similar  to  that  reached  by  Krauss  and  Stevens 
(Reference  78)  based  on  their  MCSCF  analysis  of  the  U0+  system.  Since  the 
excited  electronic  states  of  U0+  lying  in  the  region  of  strong  solar  flux  (.4 
-  .7  y)  exhibit  shifted  equilibrium  internuclear  separation  from  that  of  the 


ground  4I  state,  we  predict  efficient 
photons  for  this  system. 
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4.3  U02+. 

An  analysis  of  U0^+  has  been  carried  out  in  symmetry.  McGlynn  and 

Smith  (Reference  82)  have  proposed  a  ^  ground  state  for  this  ion  based  on 
simple  molecular  orbital  arguments.  Their  set  of  one-electron  orbital  ener¬ 
gies  is  based  on  a  maximum  overlap  criterion  that  is  empirical  in  character. 
More  modern  calculations  of  the  actinide  series  atoms  suggest  that  a 
U  5[0  ^  J  2<3r  U  3[0  structure  should  be  the  most  stable  configuration.  In 

terms  of  MO's,  the  lowest  predicted  electronic  state  (D^)  would  be: 

lag2  lau2  2ag2  liru^  2au2  3ag2  2ttu^  Hg^  3au2  |  14>u  :  2l>u  (98) 

An  extensive  series  of  calculations  of  the  possible  low-lying  symmetries  of 
U0/,+  now  definitely  establish  this  to  be  the  ground  state,  but  also  predict  a 
very  low-lying  ^A^  state. 

A  series  of  calculations  of  the  low-lying  electronic  states  of  U+^,  U+1+ 
and  U+J  was  first  undertaken  to  determine  the  approximate  location  of  the 
electronic  states  of  the  central  ion  which  corresponded  to  f  +  d  transitions. 
These  results  are  given  in  Table  3  which  indicate  that  f  d  promotion  in 
these  ions  lies  at  327,  167  and  108  nm,  respectively,  for  U+J,  U+4  and  U+i. 
Thus  strong  absorption  (dipole-allowed),  corresponding  to  central  ion  promo¬ 
tion,  is  predicted  to  occur  only  for  wavelengths  shorter  than  ~  300  nm  for 
UO^ + ,  since  the  charge  state  of  the  uranium  atom  must  be  >  +  3.  A  spin  orbit 
analysis  of  these  calculated  levels  yields  the  splittings  in  J-J  coupling. 

The  ^5/2u  state  is  the  predicted  ground  state  multiplet. 

In  order  to  establish  that  f  +  d  transitions  in  U02+  lie  at  these  short 
wavelengths,  a  series  of  molecular  calculations  of  the  low-lying  electronic 

states  was  undertaken.  The  promotion  of  an  electron  out  of  the  H  MO  of 

u 

U02+,  which  is  nearly  purely  uranium  atom  in  composition,  gives  rise  to  the 

one-electron  spectra  shown  in  Table  4.  The  u  +  u  intra-f-shel 1  transitions 

are  dipole-forbidden  and  correspond  to  the  several  weak  absorption  bands 

observed  in  solutions  of  U02  +  (Reference  83).  The  lowest  dipole-allowed  f  +  d 

transitions  on  the  central  uranium  atom  corresponds  to  the  H  ( 5f  )  +  ^A 

u  u  g 

(6d^)  excitation  at  A  ~  300  nm.  This  transition  lies  beyond  the  region  of 
efficient  solar  pumping  of  U02  +  .  Hay  (Reference  84)  reports  a  similar  transi¬ 
tion  in  his  effective  potential  studies  of  U02+,  but  at  a  somewhat  shorter 
wavelength  (~  250  nm) . 
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A  separate  molecular  spectra  for  UO  2+  arises  from  charge  transfer  states 
formed  by  promotion  of  an  electron  from  (mainly)  ligand  MO's  to  a  central 
uranium  atom  MO.  Preliminary  calculations  (Reference  76)  indicated  that 
several  of  these  states  lie  in  the  300-400  nm  region,  and  would  therefore  be 
efficiently  solar  pumped.  A  careful  re-analysis  of  the  ground  2$  state 
indicated  a  false  convergence  onto  local  minima  in  the  orbital  parameter 
space.  This  produced  a  ground  state  potential  curve  with  an  internuclear 
separation  that  was  too  large,  yielding  excitation  energies  that  were  corres¬ 
pondingly  too  small.  Careful  convergence  studies  and  the  employment  of  a  new 
integration  scheme  now  gives  an  equilibrium  U-0  bond  length  in  UO  +  of  1.73  A. 
This  value,  coupled  with  a  fundamental  vibrational  frequency  of  880  cm-1, 
is  now  in  good  agreement  with  Badger's  rule  (Reference  85)  correlating  bond 
lengths  and  force  constants  for  actinide  salts  and  oxo-ions: 

R(U-0)  =  1.08  k  _1/3  +  1.17  (99) 

Using  these  more  reliable  integration  methods,  the  low-lying  charge 
transfer  states  of  U02+  were  re-examined;  their  term  levels  are  given  in  Table 
5.  These  charge  transfer  states,  some  of  which  would  exhibit  very  strong 
absorption,  begin  at  ~  47,000  cm"1  [4n  ],  with  the  first  strong  dipole  allowed 
transitions  corresponding  to  the  excitation  1  Tig  4  1  <f<u  +  1  ng  3  l<(iu2.  These 
transitions  lie  at  53,000  -*■  57,000  cm"1,  which  correspond  to  absorption  wave¬ 
lengths  in  the  VUV,  well  outside  of  the  region  for  efficient  solar  pumping.  A 
composite  potential  energy  curve  for  U02+,  including  both  central  uranium  atom 
excitation  states  and  the  low-lying  ligand  to  central  atom  charge  transfer 
states,  is  given  in  Figure  1.  Our  calculated  assignments  of  the  f  -*■  f ,  f  ■+  d 
and  charge  transfer  excitations  among  the  low-lying  states  are  given  in  Table 


This  oxo-ion  exhibits  two  different  characteristic  equilibrium  U-0 
separations.  They  correspond  to  U+ 5  (O"2)^  configurations  with  a  short  bond 
and  to  U+4  (0"  ^  •*)  2  charge  transfer  states  which  exhibit  a  weakened  U-0  bond 
and  correspondingly  longer  equilibrium  bond  length.  The  spin  orbit  analysis 
for  U02+  is  presently  being  carried  out  but  this  analysis  will  not  affect  our 
conclusion  that  solar  pumping  of  electronically  excited  states  of  U02+  has  a 
low  efficiency. 

An  analysis  of  the  LWIR  emission  from  U02+  has  been  carried  out  based  on 
the  ground  2<j^^  electronlc  state.  These  calculations  should  also  be  repre- 


sentative  of  the  LWIR  emission  from  other  low-lying  electronic  states  since 

they  all  exhibit  similar  ionicities.  The  spectroscopic  data  for  the  . 

5/  2u 

state,  derived  from  our  calcuated  ground  state  potential  curve,  are  shown  in 
Table  7.  We  find  a  small  anharmonicity  in  the  IR  active,  asymmetric  vibra¬ 
tion,  similar  to  that  found  in  the  highly  ionic  U0+  and  UO  molecules.  Our 
calculated  f-numbers,  including  fundamentals  and  overtones  are  given  in  Table 
8  for  the  lowest  9  vibrational  levels.  We  have  also  included  the  emission 
wavelengths  for  each  transition.  For  U02+,  our  calculated  f-number  for  the 
1-0  transition  is  1.2  x  10-1+  at  A  =  11.43  p.  This  can  be  compared  with  a  cal¬ 
culated  f-number  of  5.2  x  10“ 5  at  A  =  11.3  p  for  U0+.  Relatively  strong  LWIR 
is  therefore  predicted  for  U02+. 

4.4  U02. 

An  extensive  set  of  calculations  for  the  ground  state  of  UO 2  was  carried 
out  in  D  symmetry.  The  ground  electronic  state  has  the  dominant  molecular 
orbital  configuration: 

*  E  +:  [log2  1  ou2  2og2  1  tiu^  2ou2  3og2  2  uu^  1  3  cw2  l:j>u2]  (100) 

o 

We  find  an  equilibrium  internuclear  separation  of  1.81  A,  a  value  somewhat 
larger  than  that  corresponding  to  the  various  uranyl  oxo-ions.  At  present  it 
is  uncertain  whether  this  represents  a  deficiency  in  the  calculations  or 
whether  significant  back-bonding  of  electron  charge  to  the  central  uranium 
atom  is  occurring  in  the  gas  phase  U02  species.  This  would  result  in  a 
lengthening  of  the  U-0  bond  and  a  corresponding  decrease  in  the  U-0  bond 
strength.  The  first  excited  electronic  state  of  U02  corresponds  to  a  1  <j*i  +■ 

3ttu  electron  promotion.  This  state  lies  at  ~  1.4  eV  but  we  predict  very  weak 
solar  pumping  since  the  transition  is  not  dipole  allowed.  The  first  stronglv 
allowed  transition  corresponds  to  a  1 4>u  -*■  1  6g  electron  promotion  at  ~  3.4  eV. 
This  transition  lies  at  ~  360  nm  putting  it  out  of  the  region  of  efficient 
solar  pumping. 

Calculations  of  the  electronic  structure  and  LWIR  characterist ics  of 
neutral  U02  have  been  undertaken.  Current  estimates  of  the  LWIR  properties  of 
U02  are  based  on  U02+  which  has  an  IR  active  v3  mode.  We  find  neutral  U02  to 
be  a  linear  molecule  with  geometry  similar  to  that  found  for  U02+: 
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v3  (cm"1)  880.0 

x3  (cm-1)  2.7 

r  (A)  1.73 

The  characteristic  f  Q  j  —  number  of  UO  2  (asymmetric  stretch  mode)  is  computed 
to  be  fgj  =  1.1  x  10-Lt  at  12.9  p.  This  can  be  compared  with  our  calculated 
values  for  UO  2+  ( f  ^  j  =  1*2  x  10- 4  at  X  =  11.43  g)  and  U0+  (fQ1=  3.2  x  1  0  5  at 
X  =  11.3  g)  .  The  LW1R  properties  of  UO  2  are  very  close  to  the  single  ionized 
species  U02+. 

4.5  Ti0+ . 

A  preliminary  thermodynamic  analysis  of  the  titanium/oxygen  system  lias 
been  carried  out.  The  results  are  given  in  Fig.  2  where  we  indicate  the  rela¬ 
tive  energetics  of  Ti/02  and  the  first  and  second  ionization  levels  of  this 
system.  Several  very  interesting  results  have  been  obtained.  First,  we 
confirm  that  the  associative  ionization  reaction: 

Ti  +  0  •*  TiO+  +  e  (10  1  ) 

is  exothermic  by  0.47  eV  and  would  be  an  earlv  time  (or  high  altitude 
chemistry)  conversion  route  for  Ti  +  TiO+.  TiO+  is  stable  against  further 
oxidation  by  0  atoms  hut  is  converted  to  Ti02+  by  the  reaction: 

TiO+  +  02  +  Ti02+  +  0  (102) 

which  is  exothermic  by  ~  1.6  eV.  The  highly  stable  structure  of  Ti02+ 
indicates  that  the  associative  ionization  reaction: 

Ti  +  0  2  +  Ti02+  +  e  (103) 

is  exothermic  by  ~  2.1  eV.  Thus  the  chemistry  of  the  titanium/oxygen  system 
parallels  that  of  the  uranium/oxygen  system. 

An  analysis  of  the  electronic  structure  and  spectroscopy  of  the  Ti0+  ion 
has  been  carried  out.  The  calculated  spectroscopic  data  are  given  in  Table  8. 
We  find  that  the  ground  state  of  TiO+  has  2^  symmetry.  However,  the  2ja-  statP 
lies  <  0.5  eV  in  excitation  and  several  other  low-lying  electronic  states  have 
been  found.  A  surprising  result  of  our  studies  is  that  the  dipole  moment 
function  for  the  ground  X  state  of  TiO+  is  a  weak  function  of  the  inter- 
nuclear  separation  near  Rp .  In  addition,  we  find  that  the  absolute  value  of 
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the  dLpole  moment  is  ~  7  debyes,  a  value  smaller  than  that  typically  found  for 
other  metal  oxides.  This  suggests  that  the  bonding  in  TiO+  (and  TiO)  is  more 
covalent  than  in  the  nearly  completely  filled  d-shell  metals.  We  estimate  an 
ionicity  of  approximately  73%.  Our  calculated  oscillator  strengths  for  LWIR 
in  TiO+  yields  fjg  =  4.4  x  10“®,  a  value  much  smaller  than  that  found  for  U0+ 
(5.2  x  10~5)  and  for  other  similar  metal  oxides.  The  fundamental  band  is  at 
10.3  y ,  with  a  very  weak  first  overtone  at  5.2  y.  The  oscillator  strengths 
for  the  first  six  vibrational  levels  are  given  in  Table  10.  Further  studies 
of  electronically  excited  states  of  TiO+  are  indicated  to  complete  our  LWIR 
analysis  of  this  system.  In  addition,  the  stability,  structure  and  radiative 
properties  of  Ti02+  will  be  the  subject  of  future  research. 

4.6  e  +  U+,  e  +  U++  DIELECTRONIC  RECOMBINATION. 

In  addition  to  our  calculations  of  the  structure  and  radiative  properties 
of  metal/oxygen  species,  an  analysis  of  the  importance  of  dielectronic  recom¬ 
bination  procesc.1'-  has  been  initiated.  Dielectronic  recombination  (DR)  is  the 
process  by  which  electron  capture  from  the  continuum  to  a  bound  state  is 
facilitated  by  collisional  excitation  of  a  previously  bound  target  electron. 
The  process  can  be  represented  by: 


[u+q_1] 


'[u+<*-1]  +  hv 


(104) 


(105) 


A  rt 

where  [u  ^  ]  is  a  doubly  excited  state  of  the  system  at  ionization  level 

(+q-l)  and  [u+C*  ^  ]  is  a  stabilized  excited  state  lying  below  tne  U  ^ 
ionization  threshold.  The  process  is  shown  schematically  in  Figure  3.  In  the 
case  of  recombination  unperturbed  by  external  fields  or  the  plasma  density, 
the  rate  coefficient  can  be  calculated  as: 


i ,  j  >k 


( e)  = 


^ii(e)  Ajk 

l  ^n'u)  +  V 


-Gj/kT 


where  //..(g)  is  the  rate  of  auto  ion izat ion  from  state  Ij^  into  the  continuum 
1  i 

state  |i>,  A  is  the  radiative  rate  from  state  jj>  into  the  final  bound  state 
J  k 

; k>  and  the  sum  is  over  all  possible  open  channels. 
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An  analysis  of  the  low-lying  electronic  states  of  U° ,  U+1  and  U+2  has 
been  carried  out  to  establish  a  data  base  for  examining  the  role  of  DR  in 
uranium.  The  calculated  energy  spectra  for  U^,  U+1  and  U+2  are  given  in 
Tables  11-13,  respectively.  The  states  examined  include  all  one-electron 
promotions  of  U  ^ ,  followed  by  one  electron  capture  final  configurations  which 
correspond  either  to  high-lying  Rydberg  states  or  bound  valence  excited  states 
of  the  U  ^  ^  system.  Tables  14  and  15  summarize  the  low-lying  states  which 
are  important  for  DR  of  e  +  U+  and  e  +  U++,  respectively.  There  are  two  low- 
lying  configurations  for  U++  which  have  the  J-J  coupled  representations, 

[5f  ^/2^s  ]  anc*  [5f  5/2^3/2  ] '  These  two  configurations  are  separated  by  only 

4117  cm- *  and  both  can  be  formed  by  successively  removing  two  electrons  from 

o  3  2  ++ 

the  ground  state  neutral  U  configuration,  [5f,-/27s  3  /  2  ]  ‘  ^W°  ot*ler  ^ 

representations  are  also  possible,  [  5  f  ^27s^]  anc*  [  5f^ 2  1  s  6d  3/2  ]  >  but  these 

states  lie  much  higher  > n  energy  and  would  have  negligible  population,  even  in 

a  10,000  K  plasma. 

In  the  case  of  e  +  U+  DR  (Table  14),  at  least  two  channels  are  accessible 
for  low-energy  electron  attachment.  Since  the  Rydberg  structure  lies  £  2000 
cm-1  from  the  U+  ionization  limit,  several  more  states  are  open  DR  channels 
for  electrons  with  Te  =  2000-5000  K. 

The  low-lying  excited  DR  electronic  states  of  U++  are  shown  in  Table  15, 
relative  to  both  the  e  +  U  [-^5/2^sl  and  e  +  U++ [ 5f ^ 26d ^ 2 ]  reactant 
channels.  These  data  indicate  that  the  first  open  attachment  channel  for 

*  1  •  .  ++  r  3  -I 

dielectronic  recombination  e  +  U  5fc,_7s  occurs  at  T  ~  4000  K.  However. 

L  3/2  J  e 

the  cross-section  for  capture  into  this  channel  is  very  low  since  it  involves 
a  ^3/2 — ^~^3/2  Promot  l°n  ich  has  a  low  transition  probability.  The  first 
reaction  channel  for  DR  of  e  +  U++  with  a  significant  cross-section  results 
from  the  allowed  7s — <»-6d^2  transition.  This  channel  opens  for  electron 
collision  energies,  T^  >  5000  K,  where  electron  capture  into  a  high  lying 
Rydberg  orbital  can  occur.  However,  the  probability  of  subsequent  autoioniza¬ 
tion  i.s  now  competitive  with  radiative  stabilization  and  the  dielectronic 
recombination  rate  coefficient  is  of  the  order  of  10-12  cm^/sec  or  less. 


si 


54 


We  conclude  that  dielectronic  recombination  of  e  +  U4"*'  is  a  low  probabil¬ 
ity  process  and  not  rate  competitive  with  possible  ion-neutral  charge  exchange 

reactions.  In  contrast,  dielectronic  recombination  of  e  +  U+  has  a  much 

~  3 

higher  probability  and  may  approach  gas  kinetic  rates  for  >  10  K.  Detailed 
rate  calculations  of  this  process  are  in  progress. 


4.7  N  +  02. 

As  a  final  task  under  this  Contract,  a  study  of  the  atmospheric  reactions 
0  +  N2  and  N  +  02  which  yield  vibrat ional ly  excited  NO  molecules  has  been 
initiated.  The  sources  of  NO  emission  in  the  IR  are  not  fully  understood  and 
the  detailed  kinetic  branching  of  these  reactions  is  being  studied.  In 
particular,  the  role  of  the  n[^pJ  +  02[x3£g-]  reaction  in  the  production  of 
vibrat ional ly  excited  NO  molecules  is  being  examined.  The  current  NORSE  code 
assumes  that  all  of  the  exothermic ity  of  this  reaction  goes  into  vibrational 
excitation  of  the  NO  product  molecule,  with  all  levels  being  equally  populated 
up  to  v  =  26.  In  Table  16,  we  show  the  correlations  connecting  the  N  +  02 
reactants  to  NO  +  0  product  states.  A  schematic  diagram  of  these  correlations 
for  Cg  symmetry  is  shown  in  Figures  4  and  5  for  doublet  and  quartet  states 
respectively.  This  diagram  is  a  somewhat  more  complete  version  of  that  given 
by  Donovan  and  Husain  (Reference  86). 

The  adiabatic  reaction  surfaces  arising  from  N [  *"P J  +  are  the 
following: 


n*[2p]  +  oJx3i:' 


*  r 
N02[ 

‘A'x] 

- ■•-NO* [ a4n  ] 

L  g 1 
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o[3p] 
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.212  eV 

(107) 

no2[ 

4A'vi J 

k  r 

- -NO  [a  g 

+ 

0  p] 

+ 

.212 

(108) 

k  1 

N02[ 

^A"viii J 

- ►NO  (X^Jlj 

+ 

0*[1D ! 

+ 

3.010 

(109) 

no2[ 

4  All--  I 

m  MO*  \  n  4 17  i 

0  [  3P  ] 

4. 

.212 

(110) 

A  V  J 

►  NU  [a  n  j 

T 

no2L 

3A"ix  ] 

— -no  [x^nj 

+ 

o*iisj 

+ 

.788 

(111) 

^  k  r 

no2L 

4A"vi  j 

- —NO  la4ngj 

+ 

o[3p] 

+ 

.212 

(112) 

Only  the  third  reaction  [yielding  N0(X31I)  +  0  Od)]  is  exothermic  enough  to 
yield  NO  with  a  significant  degree  of  vibrational  excitation,  provided  that 
this  reaction  proceeds  along  an  adiabatic  pathway.  However,  there  is 
spectroscopic  evidence  that  curve-c ross ing  and  non-ad iabat ic  behavior  does 
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occur  for  this  system.  The  ground  A.  state  of  NO  must  ad  iabat  ica 1 ly 

-k  2  ^  3  —  ^ 

correlate  to  the  reactants;  N  [  D]  +  E^  ],  whereas  the  first  excited  2B^ 

state  correlates  to  ground  state  reactants,  n[4s]  +  0„[x3j;  ].  In  addition, 

■  8  *r 

Gilmore  (Reference  87)  has  pointed  out  that  NO  quenching  of  both  0  [_ 1 S J  and 
0  [  1D  ]  proceeds  very  rapidly,  further  indicating  curve-crossings  in  the 
doublet  surfaces. 

In  order  to  clarify  the  branching  kinetics  of  the  N*[2P]  +  j 

reaction,  we  have  initiated  an  extensive  series  of  calculations  of  the  poten¬ 
tial  energy  reaction  surfaces  for  N02  in  the  following  symmetries:  2A' ,  2A" , 
4A' ,  4A".  Preliminary  Cl  calculations  were  carried  out  in  Cs  symmetry 
assuming  a  frozen  core  of  (la^2)  ( 1 b^ 2 )  (2a^2)  (3a^2)  (2b^2)  (4a^2),  which  is 
the  Hartree  -  Fock  representation  of  the  inner  Is  and  2s  electrons.  Unfor¬ 
tunately,  the  characteristic  charge  density  corresponding  to  this  core  repre¬ 
sentation  differs  significantly  for  different  spin  and  spatial  couplings  of  N 
+  02-  The  characteristic  charge  of  the  (4aj2)  pair,  for  example,  differs  in 
calculations  in  2A'  symmetry  as  compared  with  2A"  or  4A'  symmetry.  Thus  only 
a  frozen  inner  shell  (Is  electrons)  appears  to  be  possible  for  this  system. 
This  smaller  frozen  core  representation  results  in  a  great  increase  in  the  Cl 
size,  from  2-3000  to  20-30000.  The  calculations  must  therefore  be  carried  out 
on  a  CRAY  computer  and  the  necessary  code  modifications  for  running  this 
problem  on  the  DNA  access  machines  at  LANL  are  in  progress.  Our  calculations 
will  be  carried  out  in  Cg  symmetry  since  we  need  to  simultaneously  examine 
the  N  +  02  reaction  path  on  twelve  surfaces  in  2A'  and  2A"  symmetry  and  on  six 
surfaces  in  4A'  and  4A"  symmetry. 


SECTION  5 


RECOMMENDATIONS 

Previous  studies  of  the  electronic  structure  and  energy  levels  of  the 
uranium/oxygen  system  have  been  undertaken  by  UTRC  and  reported  in  DNA-TR-82- 
159  (Reference  75),  DNA-TR-85- 1 56  (Reference  76)  and  the  present  technical 
report  for  Contract  DNA001 -85-C-0 1 20 .  In  these  studies  the  focus  was  on 
estimating  the  LWIR  radiation  characteristics  and  the  efficiency  of  solar 
pumping  of  the  allowed  electronic  transitions  in  this  system  which  fall  in  the 
visible  wavelength  region.  An  examination  of  the  thermodynamics  and  kinetics 
of  the  uranium/oxygen  system  shows  that  U+  and  U0+  are  terminal  ions,  and 
therefore  potentially  important  radiators  at  high  altitudes,  whereas  UC>2+  is 
more  important  at  low  altitudes  where  there  is  significant  02  concentration. 
Our  preliminary  studies  also  indicated  that  TiO+  is  not  an  important  LWIR 
radiator  but  that  Ti02+  may  be  an  important  ion  at  low  altitudes.  Further 
studies  of  the  uranium/oxygen  system  are  outlined  below.  In  addition, 
suggested  studies  of  dielectronic  recombination  and  kinetic  branching  in 
atmospheric  reactions  producing  radiating  NO  molecules  are  described. 

5.1  URANIUM/ OXYGEN. 

The  calculations  of  the  electronic  structure  of  the  uranium  oxide  ions 
have  been  performed  with  neglect  of  the  spin-orbit  interaction  term, 
optimizing  the  relativistic  density  functional  SCF  wave funct  ions  for  the 
particular  ion.  Spin-orbit  interaction  is  then  introduced  as  a  perturbation 
through  conf  igurat  ional  mixing.  This  same  approach  has  been  used  by  Krauss 
and  Stevens  (Reference  78),  with  a  full  ab  initio  treatment  of  the  valence 
electrons,  using  {.-dependent  relativistic  effective  core  potentials.  Either 
approach  allows  for  polarization  and  charge  transfer  from  the  ligand  oxygen 
atoms  to  the  central  uranium  ion. 

A  different  approach  for  rare  earth  metal  oxides  has  been  described  by 
Dul  ick  (Reference  81)  which  is  based  on  a  crystal  field  model.  Application  of 
this  method  to  the  UO  and  U0+  ground  states  (Reference  89)  yields  a  signifi¬ 
cantly  different  energy  ordering  for  the  ft  coupled  states  than  is  found  in  the 
full  molecular  treatments  of  the  systems  described  above.  The  application  of 
the  crystal  field  model  to  lanthanide  monoxides,  however,  appears  to  yield  the 
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known  experimental  electronic  assignments  (Reference  90).  An  examination  of 
the  composition  of  the  predicted  ft  states  for  U0+,  based  on  this  crystal  field 
model,  indicates  a  significant  mixing  of  electronic  states  of  the  uranium  ion 
of  small  £-value.  These  effects  are  not  found  in  the  molecular  treatments, 
whose  compositions  are  dominated  by  the  tightly  coupled  5f  states.  Further 
studies  of  these  differences  in  approach  are  indicated;  however,  since  all 
methods  are  either  based  on  or  predict  a  large  ionic  character  in  the  UO+  and 
UCK+  ions,  our  conclusion  of  strong  LWIR  emission  should  not  be  altered.  The 
crystal  field  model  has  not  been  utilized  for  describing  electronic  transi¬ 
tions  or  the  absorption/emission  characteristics  in  either  the  IR  or  visible 
wavelength  regions. 

Our  suggested  studies  include  the  further  examination  of  low-lying 
electronic  states  in  U0+  and  U02+.  These  studies  will  primarily  cover  the  IR, 
visible  and  UV  excitation  wavelengths  (<  ~  5  eV)  in  order  to  more  accurately 
define  the  radiation  emission  strengths.  A  new  density  functional  computer 
code,  that  is  being  developed  in  collaboration  with  Professor  F.  E.  Harris  (U. 
of  Utah),  will  be  utilized  to  carry  out  these  suggested  studies.  This  new 
code  yields  a  more  accurate  representation  of  the  electron  density  in  the 
valence  region  between  atomic  centers.  This  region  has  previously  been 
described  by  a  constant  (muffin-tin)  potential  which  can  lead  to  inaccuracies 
in  describing  polarization  effects.  We  feel  that  this  new  code  would  also 
permit  a  more  critical  evaluation  of  the  crystal  field  model,  as  applied  to 
the  uranium/oxygen  system. 

Our  increase  in  the  data  base  for  the  uranium/oxygen  system  will  include 
radiation  f-numbers,  characteristic  absorption/emission  wavelengths  and 
calculated  density  of  states  for  U  n,  UO  n,  UO^  0  (n  =  0,1,2).  In  a  smaller 
effort,  the  radiation  properties  of  T iO+  and  TiC>2  +  should  be  more  accurately 
de  f  ined . 

5.2  DIF.LECTRONIC  RECOMBINATION. 

In  addition  to  our  suggested  calculations  of  the  structure  and  radiative 
properties  of  uranium/oxygen  species,  further  analysis  of  the  importance  of 
dielectronic  recombination  processes  should  be  carried  out.  A  preliminary 
study  of  the  important  dissociative  recombination  states  of  U°,  U+  and  U++  has 
been  described  above  under  Sect  ion  4.  Several  highly  excited  configurations 
for  11°  have  been  ident  ified  which  indicate  that  e  +  U+  may  be  an  important 


charge  neutralization  process,  but  a  more  extensive  study  is  required  to 
determine  the  density  of  states  accessible  to  dielectronic  recombination  under 
conditions  up  to  perhaps  5000°K. 

A  corresponding  DR  process  is  possible  for  e  +  U0+ .  At  the  present  time, 
only  the  optically  connected  excited  states  of  U0+  have  been  studied.  The 
ground  molecular  orbital  configuration  for  U0+  is: 

4I [lo2  2  a2  3a2  4o2  lit4  2n4  !  3  tt  1  <j>  1  6  ]  (115) 

where  thr  last  three  molecular  orbitals  are  nearly  pure  f-orbitals  centered 
on  uranium.  Excitation  processes  that  need  to  be  examined,  all  in  a  molecular 
framework,  are:  f2  ■*  d2,  f2  -*■  s2  ,  f2  ■+  sd  f2  ->•  p2,  f2  ->■  sp  and  f2  -*•  pd . 
Additional  excitations  from  the  core  2 it  level  will  be  examined  but  it  is 
expected  that  these  levels  are  too  tightly  bound  for  the  electron  energy  range 
of  interest. 

These  studies  of  dielectronic  recombination  contribute  to  our  understand¬ 
ing  of  the  kinetics  of  charge  neutralization  following  nuclear  bursts.  In 
particular,  the  role  of  dielectronic  recombination  as  a  mechanism  for  removal 
of  electrons  is  uncertain  at  present  and,  if  found  to  be  important,  must  be 
included  in  the  chemistry  kinetics  in  the  NORSE  code. 

5.3  NO/NO+  FORMATION. 

Under  this  Contract,  DNAOO 1 -85-C-0 1 20 ,  a  study  of  the  atmospheric 
reactions  0  +  and  N  +  0^  which  yield  v  ibr at  iona 1 1 y  excited  NO  molecules  was 
initiated.  This  study  was  in  response  to  the  requirement  defined  by  Gilmore 
(Reference  87)  for  a  more  accurate  analysis  of  the  n[2d]  +  0 and  n[2p]  +  0^ 
reactions  which  are  modeled  in  the  NORSE  calculation  of  NO  chemiluminescence. 
Since  the  n[2p]  +  0 2  react  ion  can  yield  n[2d]  +  0^  as  products,  as  well  as 
vibrationally  excited  NO,  there  is  a  large  uncertainty  in  the  kinetic 
branching  at  the  present  time.  If  v ibrat  iona 1 1 y  hot  NO  is  strongly  favored  as 
the  product  channel,  high  IR  emission  at  9  pm  and  4.5  pm  is  predicted.  We 
suggest  a  study  of  these  reactions,  including  all  energetically  allowed 
product  channels.  Detailed  quantum  mechanical  calculations  of  twelve  reaction 
surfaces  are  required,  necessitating  the  use  of  a  CRAY-class  computer. 

i,i  a  clos«lv  related  problem,  a  stron  radiator  in  the  IR  following  a 
nuclear  burst  is  vibrationally  excited  N0+ ,  with  emission  near  4.3  pm.  An 
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examination  of  the  ion-molecule  reaction  surface  for  0  +  Nz ,  0^  +  N,  and 

0+  +  N,  should  be  carried  out  to  establish  the  overall  product  correlations 
for  this  system.  If  required,  detailed  quantum  mechanical  studies  of  the 
pertinent  reaction  surfaces  should  be  examined.  A  similar  study  of  the 
0+  +  ■*  N0+  +  N  reaction  has  been  previously  carried  out  for  DNA  by  UTRC 

under  AFGL  Contract  FI 9628-80-C-0209  (Reference  91).  All  of  these  studies  of 
atmospheric  molecules  have  the  objective  of  increasing  the  fundamental  data 
base  of  the  reaction  kinetics  in  order  to  provide  useful  input  for  chemilumi¬ 
nescence  modeling  in  the  various  atmospheric  codes. 
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Figure  2.  Thermodynamic  stability  of  titanium/oxygen  species. 


Table  1.  Calculated  spectroscopic  data  for  IJO/UO  . 
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*SCF  and  MCSCF  results,  see  Reference  78. 

**See  Reference  79. 

***See  Reference  80. 
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Table  3.  Calculated  energy  spectra  for  u  3,  U+i+,  U+5 


1983  cm"1) 
2873  cm-1) 


2220  cm-1) 
3353  cm-1) 


2455  cm-1) 
4268  cm'M 


Configuration 

E  (cm' 

5f35/2 

0. 

5 f  2 5 / 2  5f  7/2 

6938. 

5£5/2  5f27/2 

13898. 

5f37/2 

28890. 

5f25/26d3/2 

30592. 

5f25/26d5/2 

37774. 

5E5/25f7/26d3/2 

38166. 

5f5/25f 7/26d5/2 

45342. 

5f27/26d3/2 

45754. 

5f27/26d5/2 

52938. 

5f25/2 

0 

5f5/25f7/2 

7770 

5f27/2 

15570 

5f 5/26d3/2 

60032 

3f 5/26d5/2 

68410 

5E7/26d3/2 

68912 

5E7/26d5/2 

77302 

5E5/2 

0 

5E7/2 

8592 

6d  3/ 2 

6d  r  /  -j 


92921  . 
103590. 


*ru  irviTOJ  irv  wv¥u  it.  r;  rj  »ru  #wrj iru  VJ  WV >rj »c.  V_V\.  Ww  WWAr-Tr.'VW.V.  Y.  iT./V^V-^. ^JVinlUVvVl" J^JVuV-.VJVTV 

a 

s 


Table  4.  Term  levels  for  central  atom  (one-electron)  excitation  of  U02+. 


Molecular  Orbital 


Symmetry 


Te  (cm  1) 
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32250. 
33163. 
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*Although  the  center-of-gravity  of  the  1 5u  state  lies  lower  than  1  (}u ,  spin- 
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orbit  splitting  causes  the  $5/211  state  to  become  the  ground  state,  in 
agreement  with  previous  studies  of  this  system  (Reference  83). 
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Table  10.  Calculated  oscillator  strengths  ( f  t  i, )  for  the 
vibrational  transitions  of  Ti0+. 


1 

*4.40-06 

- 

- 

- 

- 

- 

**10.26 

- 

- 

- 

- 

2 

1.87-08 

9.28-06 

- 

- 

- 

5.16 

10.35 

- 

- 

- 

- 

3 

7.46-11 

5.29-08 

1.47-05 

_ 

- 

- 

3.45 

5.20 

10.46 

- 

- 

- 

4 

2.96-11 

1.83-10 

1.03-07 

2.06-05 

- 

- 

2.60 

3.49 

5.25 

10.56 

- 

- 

5 

2.60-11 

5.63-11 

3.67-10 

1.76-07 

2.72-05 

- 

2.09 

2.63 

3.52 

5.30 

10.66 

- 

6 

2.03-13 

6.42-11 

2.56-11 

3.47-10 

2.73-07 

3.45-05 

1.75 

*Oscillator  strength 
^Wavelength,  microns 

2.11 

2.65 

3.55 

5.35 

10.76 

»J 


Table  11.  Calculated  energy  spectra  for  neutral  uranium,  U° . 


=  1832  cm 

Atomic  Configuration 

Energy  (cm- *) 

5f 3/27s26d3/2 

0. 

5f 5/27s26d 5/2 

3588. 

5e5/25E7/27s26d3/2 

6411. 

5f3/27s^7pi/2 

8691. 

5f5/25f7/27s26d5/2 

10190. 

5f2/27s26d|/2 

12815. 

5E5/25E7/27s26d3/2 

12982. 

5f5/27s27p3/2  14154. 


3  2 

5  f  ^  /  2  2  s  6d  ^  j  2 

14957. 

5E5/25E7/27s27pl/2 

15275. 

- E5/25E7/27s26d5/2 

16779. 

5f|/27s6d3/26d5/2 

18278. 

5  f  7 /2  7s2  6d  3/2 

19606. 

5E5/25E7/27s26d2/2 

19988. 

5E5/23f7/27s27p3/2 

20718. 

3E5/23E7/27s6d3/2 

21416. 

3E5/27sE>d5/2 

22000. 

5f5/25f7/27s2  7pl/2 

22191 . 

5f  7/27s26d5//2 

23343. 

2 

3E5/23E7/27sE>d3/24d5/2 

24882. 

r  c2  7  2,  .2 

f  7 / 2  7 s  6d3,/2 

27441. 

2  2 

3  E  5/2  3E7 ,/2  7s"  7 p 3 / 2 

27485. 

3  E  5  /  2  3  E  7  /  2  7  s  4d  3  / 2 

3  f  5  /  2  3  f  7 2  7  s  4d  5  /  2 

28110. 

28216. 

l;  c6d  =  1435  cm_l) 


Atomic  Configuration 

Energy  (cm 

5f 7/2  7s  2?p 1 /2 

28782. 

5E5/25E7/27s6d3/26d5/2 

31530. 

5E5/27s7pl/2 

32337. 

5f 7/27s2yp3/2 

34146. 

5E7/27s6d 3/2 

34578. 

5E5/25E7/27s6d5/2 

34857. 

3 

5f  7  ^/2^^^^3/2^^3/2 

38004. 

5E5/25E7/27s7pl/2 

38957. 

5f 7/2  7s&d 5/2 

41553. 

5E5/26d3/27pl/2 

42791. 

5E5/27s7p3/2 

44527. 

5£5/25E7/27s7pl/2 

45679. 

3  2 

5f 5/26d5/27pi/2 

47082. 

o  2 

5E5/23E7/26d3/27pl/2 

49432. 

3E5/25E7/27s7p3/2 

51368. 

5E 7/27s7p2/2 

52397. 

3E5/23E7/2'3H5/27pl/2 

53738. 

2  2 
5E5/25E7/26d3/2  7pl  ,/2 

56001 . 

3f 5/23E7/27s7p3/2 

58128. 

2  2 
3E5/23f 7/26d5/27pl/2 

60347. 

3f 7/2hd3/27pT/2 

62667. 

3  2 

3  E  7  /  2  7  s  7  p  3  /  2 

64044 . 

5  f  7  ' 2 <Sd  A  / 2  7  P  T  /  2 

67212. 

28214. 


Table  12.  Calculated  energy  spectra  for  the  uran 


ium  positive  ion,  U+ . 


Atomic  Configuration 
a  9 

5f5/27s 

5£5/23£7/27s2 

5f33/27s6d3/2 

5f|/27s6d5/2 

5£5/25£7/27s6d3/2 

5f|/27s26d3/2 

3£5/25£7/27s6d5/2 
5£5/25£7/27s 
5f 5/26d3/2 
3£5/26d3/26d5/2 
5£ |/27s6d3/2 

5f|/27s26d3/2 

3£5/27s7pl/2 

5£5/25£7/26d3/2 

5£5/26d|/2 

5f2/27s6d3/26d5/2 

5£5/25£7/26d3/26d3/2 

5£5/23£7/27s7pl/2 
5f  5/26d  3/2 

5£5/26d3/27pl/2 

5£5/27s7p3/2 

5£5/25£7/26d3/2 

3£5/27s6d5/2 

3£5/23d3/23d5/2 


(C5f  =  1840  cm"1;  U6d  =  1793  cm  1 ) 

Energy  (cm-1)  Atomic  Configuration  Energy  (cm 


0. 

5£5/26d  5/27p  1/2 

39868. 

6439. 

5£5/25£7/26d3/27pl/2 

41781. 

9589. 

5£5/25f7/27s7p5/2 

41909. 

14071. 

5£5/26d3/27p3/2 

42690. 

16151. 

5£5/26d3/26d5/2 

43397. 

19924. 

5f 5/27s27pl/2 

47345. 

20668. 

5£5/26d5/27p3/2 

47415. 

21133. 

5£5/27s6d3/27pl/2 

48506. 

21302. 

5£5/25f7/26d3/27p3/2 

49276. 

25424. 

3£5/27s6d3/27pl/2 

52645. 

25461. 

5f|/27s6d5/27p1/2 

54351. 

25582. 

3£5/27pl/2 

54553. 

26690. 

3f 5/27s  2?p3/2 

57178. 

27788. 

5f 5/27s6d3/27p3/2 

57408. 

29655. 

5£5/26d3/26d5/27pl/2 

58095. 

30689. 

5f2/26d3/27p3/2 

60825. 

31962. 

3£ 5/2  7p 1 /2  7p3/2 

63277. 

33422. 

5f?/27s6d5/,27p3/2 

63351. 

33663. 

3f 5/26d3/26d5/27p3/2 

66356. 

35050. 

3£ 5/2  7p3/2 

72295. 

2  2 
3£5/26d3/27pl/2 

77734. 

35142. 

3f 5/27s7pl/2 

78212. 

5f 5/26d3/27pl/27p3/2 

3£5/27s7pl/27p3/2 

2  2 
3£5/26d3/27p3/2 

3£5/27s7p3/2 


87191 . 
88457. 
06804. 


38568. 


90]32. 


Table  13.  Calculated  energy  spectra  for  the  uranium  positive  ion,  U++. 
(?5f  =  *932  cm  1 ;  (^6d  =  2091  cm  1 ) 


Atomic  Configuration 

5f 5/2  7s 
3 

5f5/26d3/2 

5f5/25f7/27s 

3 

5f5/26d5/2 

2 

5f5/25£7/26d3/2 

5£5/25£7/26d5/2 

5f3/27s6d3/2 

5f|/26d3/2 

3f2/27s2 

5  E 5/26d  3/26d  5/2 
5  f 5 /2  ?s6d  5 /2 

3 

5  f  5/2  7p  1  /2 

5  f  5 / 2  6d  5/2 
2 

5  f  5/2  5f  7 /2  7p  1  /2 

5  f 5/2  7p3/2 
2 

5f 5/2  5f 7/2  7p3/2 
2 

5  f  5/26d  3/2  7p  1  /2 

5C5/27s7pl/2 

2 

5  f  5  /  2  E>d  5/27pl/2 
2 

5  f  5/26d  3/2  7p3/2 
2 

5f5/27s7p3/2 

2 

5f5/26d5/27p3/2 

S f 5/2  7p  1  /2 
2 

5  f 5/2  7p 1 /2  7p  3/2 
SF2  7  2 


Energy  (cm-1) 


9345. 

10846. 


16070. 


24660. 

24722. 

27632. 

30625. 

30979. 

31198. 

36597. 

38057. 

42236. 

49092. 


57535. 


61600. 


63990. 


69400. 


74336. 


75956. 


96294. 


109336. 


122S44. 


■  «  k  «  k  »  k  *  W >  .  ■»  k^r  w 


Table  14.  Energetics  of  dielectronic  recombination  of  e  +  U+. 


e  +  U+ 


e  +  [5f|/27s2 


e  (cm  *) 


[5£ 5/27s7P1/27P3/2  ] 

[5f5/27s27p3/26d5/2  ] 
[5f 5/27s7p23/2  ] 
[5f^/27s6d3/2n£] 
[5f^/27s6d5/2n  £] 
[5f|/27s2bd3/2n£] 


e  (cm- *) 


6154. 


9589-E, 


14071-E 


19924-E_ 


■vVjflj 


-  * 


H 


vv 

vv 

Vv. 


F~FSwrfL<r.rJ.>r:~rjvr?VTJ'rJ’rJ’jvjvsvj,irr+ _r 


Table  15.  Energetics  of  dielectronic  recombination  of  e  +  U++. 


e  +  U-* 


e  +  [5fs/27s 


e  +  [5f  <^/2 6d  3/2 


e  (cm  1) 


4117. 


(cm  1 ) 


[5f 5/2 7s 7p 3/2  ] 

2  711. 

[5f 5/2^d 3/ 2  n 

4117. -En, 

[5f 5/2^d5/2n  £] 

9%5--Ent 

[5f2/27s6d3/2n  a] 

24660. -En, 

[5f  5/2bd 3/2n  «.] 

24722. -E  n£ 

[5f5/2MS/2n  *] 

5228. -En, 

[5f 3/27sbd3/2n  i] 

2°54  3.-Kn4 

[5f \i2M\/in 

20604 .  _E-n  i 

v.'-"  ~s~j- 

rv 


Table  16.  Molecular  correlation  diagram  for  N  +  0 r 


Reactants 

N(4su)  +  o2< 

N(4Su)  +  02< 

n(45u)  +  02 
n(2du)  +  o2 

N(2Du)  +  02 

N ( 2  p  u )  +  02 

n(20u)  +  02 

N(4Su)  +  0. 

N(4Su)  +  0 


w  • 

-■  .v.v 
wv 
**  %• 


:xV 


[a 


1x3 Jg> 


“V 


(x3e-) 


(blEp 


Energy  (eV) 

States 

Point 

Group 

1.401 

2 , 4 , 6  £+ 

C«v 

2,4,6B 

D  1 

C2v 

2,4,6a, 

C 

s 

HI 

2.383 


3.037 


3.785 


4.766 


4.977 


5.420 


4.  4r 
A2>  B1 

4a" , 4a' 

4  r 


4a" 


2,4^,  2,4 j|,  2,4a 
2>4B1(2),  2’4B2,  2>4A2,  2’4A1 


2  AH 


(2), 


2>V 


(3) 


1  ? ,  2  r,  2  n,  2  a,  2  4>,  2  r 

2Al(2),  2A2(3),  2B1(3),  2B2(2) 
2A'(5),  2 A" ( 5 ) 

2’4zr,  2,4n 

2,4a  2,4a  2,4b0 

a2>  2 

2  ,4A'  ,  2  ’  4A"  ( 2 ) 

2  r,  2  n,  2a 


'2v 


c«V 


'2v 


«v 


'2v 


«v 


'2v 


'2v 


»  ,J2 
2 ’ 4 ’ 6 a ’  2  > A 


J2v 


■ -SJ 


:wf 


sgss0 

;->i5 


2Aj,  2A2(2),  2 B| ,  2B2 

C2v 

.  •  -r  i 

•  >  .*» 

2A’(2),  2 A" ( 3 ) 

C 

S 

5.499 

V 

C«v 

•  ] 

% 

C2v 

V 

Cs 

5.701 

2>4'61 

2,4,6.  2 , 4 , 6„ 

C«v 

• 

\  .  N  .V 

.  'C.  - 


81 


Table  16.  Molecular  correlation  diagram  for  N  +  0 2  (Concluded). 

Point 

Reactants 

Energy  (eV) 

States 

Group 

n(4su)  +  o2(a3z*) 

5.  789 

2  >  4, 6  £- 

C«v 

2>4>6b2 

C2v 

2 , 4  >6^ti 

Cs 

N(2pu)  +  OjCa^g) 

5.958 

2  n  2  a  2  . 

n,  A,  $ 

C*v 

2A1(2),  2A2,  2B1>  2B2(2) 

C0 

2v 

2A '  ( 3 )  ,  2 A"  ( 3 ) 

C 

s 

N(2Pu)  +  02(b1 E~) 

6.612 

2r,  2n 

C«v 

2Ap  S,  2B2 

C2v 

2A'(2),  2A" 

Cs 

Point 

Products 

Energy  (eV) 

States 

Group 

no(2ii)  +  o(3pg) 

0 

2,4  e*.  2,4r,  2,4n,  2,4a 

C«v 

2 ,4A' ( 3 ) ,  2 ,4A" (3) , 

Cs 

NO(2n)  +  0( 1 D  ) 

1.967 

2  r4",  2  r,  2  n(2) ,  2  a,  2  $ 

C«v 

2 A ' ( 5  )  ,  2 A" (5) 

Cs 

NO(2n)  +  o(Lsg) 

4.189 

2  n 

C«v 

2A' ,  2A" 

Cs 

N0( a4  n)  +  0(3Pg) 

4.  765 

2,4,4 xr4,  2,4,6 r,  2,4,6n,  2 ’ 4 ’ 4 a 

c 

«V 

NO(A2£+)  +  0  ( 3Pg) 

5.450 

2 ’ 4 ’  6A ' ( 3 ) ,  2 ,4» 6A" ( 3 ) 

2,4 r,  2,4n 

eV, 

2’4A'  ,  2  ’ 4  A"  (  2  ) 

Cs 

NO(B2n)  +  0(3Pg) 

5.693 

2,4e+,  2,4r,  2,4n.  2,4a 

c 

OV 

2 ,4A 1 ( 3 ) ,  2 ’ 4A" ( 3 ) 

Cs 

NO( b4  Z~ )  +  °(3pg) 

6.0344 

2,4,6E+;  2 ,4,6  n 

c 

«v 

2,4,6A'(2),  2,4,6A" 

c 

fi 

82 
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